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MODELING  OF  DIFFUSE  PHOTOMETRIC  SIGNATURES 
OF  SATELLITES  FOR  SPACE  OBJECT  IDENTIFICATION 

Abstract 

The  diffuse  reflective  characteristics  of  four  types  of  low  earth 
orbit  satellites  were  mathematically  modeled  using  phase  functions  for 
Ideal  Lambertian  surfaces.  A  FORTRAN  computer  program  was  developed  to 
generate  simulated  signatures  and  compare  them  point  by  point  to  real 
signatures  to  obtain  a  sum  of  the  squares  of  the  residuals  (SSR),  in 
order  to  perform  pattern  recognition  and  satellite  Identification. 

Photometric  signatures  collected  by  the  Satellite  Identification 
and  Tracking  Unit  (SITU),  St.  Margarets,  Canada  wer^recelvetf^from  the 
ADC  Intelligence  Center  In  thefUORAD  Cheyenne  Mountain  Complex  In 
Colorado  Springs,  Colorado,  atane  set  of  signatures  was  used  to  validate 
the  computer  model  of  one  satellite  type,  and  the  others  were  used  to 
test  the  program's  ability  to  Identify  the  satellite.  The  tested  model 
Involved  llne-of-slght  obcuratlon  of  same  parts  by  others,  relative 
motion  of  body  parts,  and  for  phase  shadowing. 

The  program  was  able  to  correctly  Identify  the  modeled  satellite, 
as  long  as  the  phase  angle  remained  small,  generally  less  than  ninety 
degrees.  For  larger  phase  angles,  the  true  signatures  diverged  signi¬ 
ficantly  from  simulated  signatures.  In  every  case,  the  signatures 
predicted  using  the  Lambertian  assumption  were  dimmer  than  the  measured 
signatures  at  larger  phase  angles. 
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Photometric  pattern  recognition  of  satellites  using  phase  function 
models  appears  to  be  feasible,  but  satellite  models  contained  In  an 
operationally  useful  computer  program  must  be  valid  for  any  viewing 
geometry,  and  should  therefore  account  for  the  non-Lambertlan  behavior 
of  Illuminated  surfaces  where  viewed  at  large  phase  angles,  or  away 
from  normal  Incidence. 


I.  INTRODUCTION  AND  BACKGROUND 


Photometric  Space  Object  Identification 

Space  Object  Identification  (SOI)  The  secondary  missions  of  the 
North  American  Aerospace  Defense  Command  (NORAD)  are  space  track,  to 
maintain  current  orbital  elements  and  provide  reentry  predictions  for 
all  man-made  space  objects,  and  SOI,  to  determine  the  physical  and 
dynamic  characteristics  of  satellites  In  near  earth  or  deep  space  orbits 
(Ref  21).  SOI  sensors  include  the  majority  of  the  NORAD  space  track 
radars,  two  contractor  operated  wideband  coherent  radars,  and  two  sat¬ 
ellite  tracking  astronomical  observatories  which  employ  photoelectric 
photometers  to  collect  time  vs.  amplitude  plots,  or  signatures,  of  the 
sunlight  reflected  from  the  satellites.  One  facility  also  collects 
photographic  Images  of  satellites  in  low  earth  orbit  (altitude<1000km). 
This  paper  is  solely  concerned  with  the  analysis  of  visible  light  photo¬ 
metric  signatures  of  satellites  and  their  application  to  the  NORAD/ADCOM 
SOI  mission.  Table  1-1  describes  the  two  currently  operating  dedicated 
photometric  SOI  sensors,  and  some  characteristics  of  the  soon  to  be 
operational  Ground-Based  Electro-Optical  Deep  Space  Surveillance  (GEODSS) 
system  (Ref  10,  Ref  11).  Table  1-2  lists  all  of  the  SOI  sensors. 

The  Application  of  Satellite  Photometry  to  SOI  No  single  SOI  data 
type  can  provide  an  exhaustive  description  of  all  satellite  character¬ 
istics  of  Interest,  but  an  Integrated  approach,  employing  Inputs  from  a 
variety  of  sources  is  necessary  (Ref  16:6).  Some  contributions  of 


A«V^. 


photometric  signatures  to  the  overall  picture  are  unique. 

Table  1-3  lists  some  of  the  Information  content  of  photometric 
signatures  compared  to  radar  signatures  (Ref  7,  Ref  14). 

Photometric  Analysis  Capability  At  the  AD1C 

In  1977*  the  General  Electric  Company,  under  SAMSO  contract,  pub¬ 
lished  the  results  of  a  system  requirements  and  functional  description 
study  of  a  proposed  SOI  Central  Analysis  System  (SOICAS)  (Ref  19  and 
Ref  20).  A  reason  for  the  SOICAS  studies  was  recognition  that,  at  that 
time,  SOI  analytical  capabilities  In  general  were  Inadequate  to  satisfy 
either  the  tactical  time  constraints  of  the  ADCOM  space  defense  mission 
or  the  data  analysis  resolution  requirements  (Ref  20:  1-1-1).  While 
great  strides  have  been  made  since  then  In  the  analysis  and  processing 
of  radar  data,  the  situation  has  remained  essentially  unchanged  for 
photometric  analysis,  despite  the  acquisition  In  1977  of  an  operational 
photometric  analysis  software  package,  the  Photometric  Data  Analysis 
Nodule  (POAM),  developed  by  the  AVCO  Systems  Division  of  AVCO  -  Everett 
Research  Laboratories  (Ref  7:1v). 

POAM  Capabilities  According  to  the  SOI  analyst  PDAM  training  course. 
"PDAM  Is  an  Interactive  software  package  Integrated  on  the  HIS  6060  com¬ 
puter  In  the  ADIC,  which  has  as  Its  main  purpose  the  analysis  of  photo¬ 
metric  signatures  (specular  and  diffuse),  from  which  estimates  of  tar¬ 
get  configuration/shape,  size,  orientation,  surface  properties  (and) 
motion  are  obtained  (Ref  14:18).N  The  software  package  contains  three 
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Mount  Equatorial  on  an  Azimuth  4-Axis  Modified 

_ Tahle^  Hydrostatic  bearings _ Baker-Nunn 
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operating  phases  for  photometric  signature  analysis  (Ref  14:40-41). 

The  preprocessing  phase  enables  an  analyst  to  retrieve  any  photo¬ 
metric  signature  from  the  Intelligence  Data  Handling  System  (IDHS) 
signature  wraparound  file,  or  from  a  disk  containing  permanently  stored 
data.  The  signature  may  be  displayed  In  whole  or  in  part  with  the 
vertical  (stellar  magnitude)  scale  selected  by  the  analyst.  The  ana¬ 
lyst  may  edit  out  portions  of  the  signal  not  desired  for  analysis  (noisy 
parts,  stars,  gaps  In  data)  and  save  portions  for  later  specular  and 
diffuse  analysis  (Ref  14:40-41).  In  practice,  this  portion  of  the  PDAM 
program  has  proven  to  be  extremely  valuable  and  easy  to  use. 

Specular  and  Diffuse  analysis  make  up  the  second,  and  most  Important 
phase  of  a  PDAM  analysis  (Ref  14:40).  The  specular  and  diffuse  reflec¬ 
tion  components  of  a  signature  are  analyzed  separately  to  yield  very 
different  kinds  of  Information,  and  the  two  signature  components  are 
separately  stored  by  the  analyst  during  preprocessing  (Ref  14:33-34). 

The  specular  signature  component  Is  used  to  determine  the  rotational 
period  of  spinning  or  tumbling  objects,  the  surface  area  of  specularly 
reflecting  surfaces  (flat  plate  or  cylinder),  curvature  of  specularly 
reflecting  spherical  segments,  the  orientation  of  surface  normals,  and 
angular  momentum  vector  direction  (Ref  14:147  and  Ref  15).  The  diffuse 
components  can  be  used  to  obtain  an  estimate  of  an  object's  size,  shape, 
orientation,  rotational  period  and  angular  momentum  vector  direction 
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The  third  phase  of  a  PDAM  analysis  Is  the  results  summary,  which 
provides  a  detailed  tabulation  of  the  results  of  all  specular  and  dif¬ 
fuse  analysis  performed  on  a  particular  signature.  Including  a  graphics 
line  drawing  of  the  satellite  model  obtained  from  the  analysis  (Ref  14: 
40). 

Operational  Limitations  of  PDAM 

Operational  experience  with  PDAM  has  shown  It  to  be  an  effective 
analytical  tool  for  specular  analysis,  but  It  rarely  yields  SOI  Informa¬ 
tion  through  diffuse  analysis  which  could  not  be  more  easily  or  quickly 
extracted  from  radar  signatures,  at  least  for  low  earth  orbit  satellites. 
The  reason  for  this  becomes  clear  when  we  compare  the  PDAM  approach  to 
diffuse  analysis  with  real  world  SOI  analysis  requirements. 

PDAM  Diffuse  Analysis 

The  diffuse  analysis  module  allows  the  analyst  to  select  a  satellite 
model  from  a  shape  library  containing  18  simple  combinations  of  shapes, 
each  of  which  has  well  understood  diffuse  light  scattering  characteris¬ 
tics  (Ref  14:134).  Figure  1-1  depicts  the  PDAM  shape  library  (Ref  14: 
265-282).  Each  library  shape  Is  characterized  by  a  number  of  parameters 
which  specify  the  dimensions  and  reflectivities  of  features,  and  the 
orientation  of  the  body  as  a  whole  In  either  an  Inertial  or  a  body-fixed 
coordinate  frame  (Ref  7:81).  The  parameter  values  selected  by  the  analyst 
for  a  given  shape  constitute  the  parameter  vector,  P  for  the  diffuse 
satellite  model  (Ref  7:81). 


Shape 

Shape 

#  of  P-Vector 

Number 

Description 

Parameters 

Illustration 

1 

Sphere 

1 

O 
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Flat  Plate 

3 

CJ 
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Cylinder 

4 

D 
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Cone 

4 

O 
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Cyl -Pit 

5 

o~) 

6 

Cone-Pit 

5 

<□1 

7 

Cone-Cyl 

5 

8 

Cone-Plt-Cyl 

6 

Qr> 

9 

Cyllnder-Spher  Endcaps 

4 

(Tl) 

10 

Sphere-Pit 

4 

*0" 

11 

Cyl-1  Frustum 

6 

(O 

12 

Cyl-2  Frustums 

6 

13 

Cyl -Rocket  Nozzle 

6 

Sri 

14 

Cyl-Frust-Noz 

7 

o 

15 

Sphere* 

2 

© 

16 

Cyl -Flat  Side  Pit 

8 

17 

Cyl-Edge  Side  Pit 

8 

18 

Cyl -Paddle 

8 

v 

♦This  sphere  differs  from  shape  l  In  that  the  analyst  can  vary  the 
ratio  of  specular  to  diffuse  reflectivities.  Model  1  Includes  only 
diffuse,  assuming  a  perfect  lambertlan  sphere. 

Figure  1-1 

The  PDAM  Shape  Library 
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The  analyst  uses  a  three-fold  approach  to  diffuse  analysis,  begin¬ 
ning  with  choice  of  a  library  shape.  If  there  Is  no  a  priori  knowledge 
of  the  satellite  shape,  all  18  of  the  library  shapes  should  be  chosen 
In  turn  (Ref  7:80).  Once  the  P-vector  components  are  specified,  the 
program  executes  two  Iteration  loops,  or  three.  In  the  case  of  unknown 
satellite  orientation,  and  calculates  the  best  parameter  values  for  the 
chosen  shape  through  successive  comparison  of  measured  signature  data 
points  to  data  points  synthetically  calculated  from  the  satellite  model 
composite  phase  function  (Ref  7:36).  The  program  then  calculates  a 
weighted  mean  squre  error  criterion  for  the  chosen  shape.  The  analyst 
wishes  to  minimize  the  sum  of  the  squares  of  residuals  function  (SSR), 


given  by 


n  2 

ssr  =  -y- 

n  Is  the  number  of  observations, 

M(tj)  Is  the  1th  data  point  of  the  measured  signature, 
S(t^P)  Is  the  1th  synthetic  data  point  (Ref  14:132),  and 
<T  Is  given  by 


U-dJ 


where 


cL  - 


The  program  calculates  the  gradient  of  the  SSR  function  with  respect  to 
each  parame.er  and  Iteratively  solves  for  each  parameter,  (Ref  14:132): 

4=2  z[s(U)  -M  tt0][-~— ] 

An  iterative  process  Is  used  instead  of  equating  grad  (SSR)  to  zero  and 
solving  forp,  because  S(t^P)  has  a  non-linear  dependence  on p (Ref  14: 
132). 

The  new  P  vector  Is  given  by 

P«ru)  '  P 

The  function  minimization  algorithm  used  here  Is  the  Jacobson  Gradient 
Method  (Ref  7:81 -A).  Finally,  an  error  matrix  describing  error  associat¬ 
ed  with  each  of  the  analyst's  Input  parameters  and  the  total  SSR,  includ¬ 
ing  convergence  or  non-convergence  of  the  Iteration  process  Is  output  In 
a  results  summary  (Ref  14:132). 

The  analyst  repeats  the  above  process  for  a  single  library  shape, 
altering  the  input  parameters  until  a  converging  solution  with  minimum 
SSR  for  that  shape  Is  obtained,  or  until  It  Is  determined  that  the  chosen 
shape  will  not  allow  a  convergent  solution. 

The  third  step  In  PDAM  diffuse  analysis  Is  selection  of  the  shape 
which  yields  the  minimum  SSR  with  respect  to  the  real  data,  the  P-vector 
components,  being  the  ones  which  produced  the  best  data  fit  for  the 
chosen  shape  (Ref  7:81).  Figure  1-2  sumnarlzes  the  diffuse  analysis 
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model  parameter  estimation  process  (Ref  18:16  and  Ref  14:131). 

Intrinsic  Limitations  of  Diffuse  Analysis  The  diffuse  signature 
of  a  satellite  is  dependent  upon  sun-target-observer  geometry,  target 
size  and  shape,  the  reflective  properties  of  target  surface  materials, 
and  target  dynamics  (Ref  14:111).  Each  of  these  factors  imposes  limi¬ 
tations  on  the  utility  of  diffuse  data. 

The  combination  of  sun-target-observer  geometry  and  target  dynamics 
has  the  greatest  influence  on  a  diffuse  signature.  The  easiest  case  to 
analyze  is  the  rota tin^sa tell ite  in  low  earth  orbit.  From  Figure  1-3 
we  can  see  that  the  greatest  time  rate  of  change  of  sensor  aspect  angle, 
otand  phase  angle, ^  ,  occurs  for  the  low  earth  orbit,  spinning/tumbling 
case.  This  greatly  reduces  the  possibility  that  more  than  one  set  of 
model  parameters  will  produce  a  signature  which  closely  matches  the  real 
data  (Ref  14:19).  The  most  difficult  case,  on  the  other  hand,  is  the 
stable  satellite  in  geosynchronous  orbit.  Phase  angle  and  aspect  angle 
both  change  very  slowly,  and  many  redundant  solutions  are  possible  (Ref 
14:19).  In  ROAM,  many  parameter  sets  may  produce  a  converging  solution 
with  respect  to  the  signature  of  a  stable  object  (Ref  7:105).  Table 
1-4  summarizes  the  effect  of  different  geometries  on  the  difficulty  of 
achieving  a  good  solution  (Ref  14:128). 

When  neither  the  size  nor  the  reflectivity  of  a  target  are  known, 
these  two  parameters  must  be  taken  together  as  a  product  of  reflectivity 
and  projected  area,  unless  multispectral  data  are  available  to  determine 
reflectivity  (ref  7:81). 


The  diffuse  analysis  limitation  which  is  most  responsible  for  the 
previously  mentioned  limited  usefulness  of  PDAM  in  the  operational 
environment  becomes  obvious  upon  inspection  of  Figue  1-1,  Many  actual 
satellite  shapes  cannot  be  closely  approximated  by  any  member  of  the 
shape  library.  In  general,  the  simple  shape  library  approach  is  accur¬ 
ate  for  very  simple  satellites,  but  is  useful  only  for  obtaining  esti¬ 
mates  of  gross  shape  and  orientation  of  more  complex  objects  (Ref  18:78) 
unless  the  models  are  more  sophisticated.  When  the  true  satellite  con¬ 
figuration  does  not  closely  match  a  library  shape,  diffuse  analysis 
results  may  be  meaningless.  Determining  the  size,  shape  and  altitude 
of  a  satellite  from  the  photometric  time  history  of  the  brightness  of 
an  unresolved  point  source  is  further  complicated  by  noise  in  the 
observations,  the  lack  of  true  Independence  between  separate  observa¬ 
tions,  and  the  impossibility  of  fully  describing  a  complex  object 
uniquely  with  a  few  parameters  (Ref  18:59). 


— Dynamics 

Orbit 

Stable 

Tumbling  or 
Precessing 

3-Axis  or 

Earth  Oriented 

Spin 

low 

earth 

(ALT  1000km) 

large  d+/dt 
small  dm/dt 
large  dp/dt 
intermediate 
difficulty 

same  as  3-axis 
but  also  good 
for  specular 
analysis 

large  d+/dt 
large  dm/dt 
large  dp/dt 
easiest  diffuse 
case 

deep 

space 

(ALT  1000km) 

small  d4/dt 
small  d«/dt 
small  dp/dt 
long  tracks 
possible 
toughest 
diffuse  case 

same  as  3-axis 
but  also  good 
for  specular 
analysis 

small  d$/dt 
large  do/dt 
large  dp/dt 
case  of  intermediate 
difficulty.  Long 
tracks  possible 

<f>  is  phase  angle,  the  angle  between  the  sun  direction  and  line  of 
sight. 

o<  is  sensor  aspect  angle,  the  angle  between  the  satellite  longi¬ 
tudinal  axis  of  symmetry  and  the  sensor  line  of  sight 
p  is  solar  aspect  angle,  the  angle  between  the  satellite 
longitudinal  axis  of  symmetry  and  the  sun  direction 


Table  1-4 

Effects  of  Orbit  and  Satellite  Dynamics  on 
the  Difficulty  of  Diffuse  Analysis 


Real  World  SOI  Requirements  vs.  PDAM  Capabilities 


At  present,  there  Is  rarely  a  need  to  use  PDAM  to  obtain  gross 
estimates  of  the  size,  shape  and  orientation  of  unknown  objects  because 
other  data  sources  are  available  which  yield  more  definitive  results. 
The  utility  of  PDAM  In  this  role  will  probably  increase,  however,  when 
the  Eastern  hemisphere  GEODSS  sensors  become  operational.  In  the  mean¬ 
time,  PDAM  is  not  useful  as  an  identification  aid  for  known  types  of 
satellites  because  of  the  generality  of  its  shape  library. 

When  the  GEODSS  system  comes  on  line,  it  Is  likely  that  the 
opportunities  for  an  electro-optical  sensor  to  detect  newly  launched 
objects  will  increase,  and  photometric  analysis  software  capable  of 
performing  pattern  recognition  for  early  mission  Identification  could 
be  very  useful.  Such  software  might  also  be  useful  for  Identification 
of  uncorrelated  targets  (Ref  23,  Ref  10). 

The  early  research  In  satellite  photometry,  which  laid  the  ground¬ 
work  for  the  PDAM  program,  was  based  upon  a  choice  between  two  possible 
approaches  to  analysis  of  photometric  satellite  signatures  (Ref  18:60). 
The  approach  chosen  for  PDAM  was  to  consider  a  family  of  simple  shapes 
and  combinations  of  shapes,  and  to  ask  what  the  dimensions  and  reflec¬ 
tivity  of  the  selected  shape  would  have  to  be  to  produce  the  observed 
signature,  beginning  with  no  a  priori  knowledge  of  the  true  satellite 
shape  (Ref  18:60). 


The  other  approach  Is  to  consider  a  family  of  known  target  satel¬ 
lite  shapes,  and  to  predict  the  photometric  signature  which  would  be 
observed  by  a  sensor  If  the  known  shape  selected  were  an  accurate  model 
of  the  satellite,  and  to  compare  the  actual  data  to  predicted  signatures 
for  all  of  the  known  shapes.  In  order  to  select  the  best  match  (Ref  18:60). 
Such  an  approach  amounts  to  pattern  recognition,  which  may  be  loosely 
defined  as  "any  automatic  system  which  does  tasks  labeled  detection,  rec¬ 
ognition,  Identification  or  classification  (Ref  16:2)."  Gamache  and  La 
Rosa  recognized  the  value  and  feasibility  of  this  approach  In  their  sig¬ 
nature  prediction  study  of  cylindrical,  solar  cell  covered  geosynchron¬ 
ous  satellites  (Ref  13:209-211). 

Statement  of  the  Problem 

This  thesis  proposes  that  both  of  the  above  approaches  to  diffuse 
analysis  are  necessary  to  fulfill  the  MJCOR  photometric  data  analysis 
mission  at  the  AD1C,  and  seeks  to  demonstrate  that  photometric  pattern 
recognition  Is  possible  for  stable  payloads  which  yield  purely  diffuse 
signatures. 

Scope  of  the  Project 

The  goals  of  this  thesis  project  are  to: 

1.  Mathematically  model  the  diffuse  light  scattering  characteristics 
of  several  stable  foreign  payloads  with  sufficient  accuracy  that  the 
models  can  be  used  for  pattern  recognl tlon  purposes. 
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2.  Validate  and  verify  the  diffuse  models  using  open  source  and 
intelligence  information  to  obtain  accurate  physical  and  dynamic 
characteristics,  and  statistical  comparisons  with  actual  signatures 
collected  by  the  NORAD  photometric  sensors. 

3.  Write  a  FORTRAN  program  which  will  automatically  compare  a  real 
signature  to  a  series  of  synthetic  signatures  for  each  satellite  model 
,  and  select  the  model  which  produces  the  signature  most  closely 
matching  the  original. 

The  project  Is  restricted  to  diffuse  data  on  stable  objects  In 
low  earth  orbit,  because  the  stable,  diffuse  case  Is  of  the  greatest 
Intelligence  Interest,  and  because  more  signatures  and  satellite 
configurations  are  available  for  low  earth  orbit,  thus  making  valid¬ 
ation  of  the  models  easier. 

General  Approach 

1.  Collect  photometric  signatures  with  their  position  and  velocity 
vectors.  These  were  obtained  with  the  cooperation  of  the  photometric 
analyst  at  the  ADIC. 

2.  Write  an  orbit  prediction  and  sensor  look  angles  computer  prog¬ 
ram. 

3.  Develop  diffuse  satellite  scattering  models  and  Incorporate  them 
Into  the  computer  program  to  generate  synthetic  signatures. 

4.  Add  a  statistical  data  comparison  subroutine  to  the  computer 
program. 


5.  Digitize  NORAD  photometric  signatures  and  use  them  to  validate 
satellite  models. 

6.  Use  the  completed  computer  program  to  automatically  Identify 
satellite  signatures. 

Sequence  of  Presentation 

The  thesis  Is  organized  as  follows: 

1.  Introduction  and  Background,  Chapter  I 

2.  Theory  of  Satellite  Photometry,  Chapter  II 

3.  Brief  Functional  Description  of  the  Computer  Program, 

Chapter  III 

4.  Satellite  Models  and  Validation  Results,  Chapter  IV 

5.  Results  of  the  Pattern  Recognition  Experiment,  Chapter  V. 

6.  Conclusions  and  Recommendations,  Chapter  VI. 

7.  Appendices 


II.  THEORY  OF  SATELLITE  PHOTOMETRY 


Photometry  Is  the  measurement  of  the  Inradiance  of  light  emitted 
from  a  source.  Applied  to  earth  satellites  It  Is  "the  measurement  and 
Interpretation  of  the  solar  energy  reflected  from  an  orbiting  body  and 
the  time  variation  of  the  reflected  energy  (Ref  7:5). "  Satellite 
photometry  employs  the  Instrumentation  and  data  gathering  techniques 
of  astronomical  photometry,  and  suffers  from  the  same  environmental  and 
Instrumental  sources  of  measurement  error.  The  principal  difference  be¬ 
tween  the  two  Is  the  nature  of  the  objects  observed. 

The  Intensity  of  sunlight  reflected  from  a  satellite  Is  dependent 
upon  sun-satelllte-sensor  geometry,  the  reflectivities  of  satellite 
components  and  the  satellite's  motion.  Measurement  of  this  Intensity 
Is  affected  by  the  atmosphere,  optics  and  electronics  associated  with 
the  photometer  (Ref  9:2).  The  high  angular  velocity  of  satellite  against 
the  sky  background  Is  an  additional  complication.  The  astronomer  need 
be  concerned  only  with  slow  changes  In  line  of  sight  air  mass,  sky  bright 
ness  and  atmospheric  stability  during  the  period  of  observation.  A 
satellite  however,  may  cross  many  degrees  of  sky,  resulting  In  a  rapidly 
changing  background.  Table  II-l  lists  the  factors  Influencing  satellite 
signatures  (Ref  7:19-26,  Ref  9:9,  Ref  3:76,  408). 

Equivalent  Stellar  Magnitude  Satellite  photometry  employs  the 
astronomical  stellar  magnitude  scale,  originated  by  Hipparchus  over 
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TABLE  II-l  Factors  Influencing  Photometric  Signatures 
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2000  years  ago,  and  in  common  use  by  astronomers  since  the  second 
century  A.O.  (Ref  2:5).  In  1856,  N.R.  Pogson  standardized  the  scale 
by  defining  a  difference  of  5  magnitudes  as  equalling  a  difference  in 
photon  flux  of  100  times  (Ref  2:5).  Therefore,  a  difference  in  Intensi¬ 
ty  of  one  magnitude  is  the  same  as  a  flux  ratio  of  100*^=10^.  Com¬ 
paring  the  intensities  of  two  objects  of  magnitudes  and  Mg,  if 
Mg  -  Mj  =  n  magnitudes,  the  Irradiance  ratio 

-£*(10*)"  ,  ••• 

Mi  *  - 2  5  l*v  (r.) 

(Ref  2:5,25)  In  satellite  photometry  magnitude  Mg  becomes  the  exoatmos- 
pherlc  solar  magnitude  (-26.78),  Ej  is  the  observed  Irradiance  from  the 
satellite  and  Eg  Is  the  solar  irradiance  in  the  visible  bandpass.  The 
expression  for  equivalent  stellar  magnitude  of  a  satellite  is 

M=  -2.S  l«3f0  £  *s 

reflected  irradiance,  and  EQ  Is  solar  Irradiance  in  the  visible  bandpass 
(  616  w/in2  from  3800  to  7600  A)  (Ref  7:6-12).  In  the  stellar  magnitude 
scale,  greater  brightness  Is  represented  by  smaller  numbers.  The  bright¬ 
est  star  in  the  sky  other  than  the  sun  (slrius)  has  a  magnitude  of  -1.58. 
The  faintest  stars  visible  to  the  naked  eye  have  a  magnitude  of  about 
+  6.0  (Ref  7:4). 


Observed  Irradiance  The  total  irradiance  of  an  object  seen  through 
the  atmosphere  in  a  bandpass  \  to  \ is  given  by  (Ref  2:26), 

ru 

E-  I  'UtM  T.lMTf  dX 

•  X 

where  T*(A)  is  atmospheric  transmission, 

T,(A)  is  transmission  of  the  optics 
T^(A)  is  transmission  of  filters,  and 
E(A)  is  exoatmospheric  flux. 

Sometimes,  all  of  the  quantities  in  the  integral  having  to  do  with 
optics  and  the  detector  are  combined  into  an  overall  detector  response 
function, ft (A),  leading  to  the  expression  (Ref  9:4), 

r' 

E«  I  T*(A)  R(A)  E(A)  dA 


This  project  deals  with  data  which  have  been  corrected  for  atmos¬ 
pheric  extinction  and  background  radiation  at  the  sensor,  so  satellite 
signature  predictions  are  for  the  exoatmospheric  case. 

The  final  general  expression  for  observed  irradiance  becomes  (Ref 


7:12), 


A 

. 


UX)  fcMd* 


Diffuse  Reflected  Irradiance  This  section  presents  the  basic 
principles  of  diffuse  reflection  and  introduces  the  diffuse  phase  func¬ 
tions  of  the  simple  shapec  which  are  combined  to  form  the  satellite 
models. 

Diffuse  reflection  occurs  when  a  surface  has  irregularities  which 
are  large  with  respect  to  the  wavelength  of  incident  radiation  and  the 
radiation  is  reflected  in  all  directions,  or  isotropically  (Ref  4:65). 
Most  surfaces  exhibit  both  diffuse  and  specular  reflection  to  some 
degree.  Figure  II-l  illustrates  specular  and  diffuse  reflection  prop¬ 
erties. 

Lambert  defined  a  perfect  diffuse  reflector  as  one  which  has 
constant  radiance,  L  (w/M  -sr),  regardless  of  the  angle  of  reflection 
from  the  surface  normal  (Ref  3:95).  In  other  words,  the  radiance  of 
the  surface  as  seen  by  a  sensor  is  not  dependent  upon  the  viewing  angle, 
if  the  surface  fills  the  entire  field  of  view  (Ref  3:531).  Since  sat¬ 
ellites  are  at  great  distances  from  earth  based  sensors  they  do  not 
fill  the  entire  field  of  view  and  may  be  considered  point  souces  for 
photometry  (Ref  3:535-537). 

Point  source  irradiance  Is  governed  by  the  inverse  square  and  cos¬ 
ine  laws  of  irradiance.  If  P  is  a  point  source  of  radiant  intenstlty  I 
and  distance  r  from  some  surface  element  dA,  and  the  normal  to  dA  Is  an 
angle  6  from  the  direction  of  P,  then  the  radiant  flux  incident  over 
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Figure  II-l  Specular  and  Diffuse  Reflection 


25 


cl#  -  I  c Ul 


where  dJl  is  the  solid  angle  subtended  by  dA  from  P  (Ref  3:93).  By  the 
definition  of  solid  angle. 


if  the  normal  to  dA  lies  along  a  radius  r,  or 

dil=  dA_c^s£ 

rl 

if  the  normal  is  an  angle  0  from  the  direction  of  P  (Ref  3:94).  By 
definition,  irradiance,  E  is  given  by  (Ref  3:92). 

3  A 

so  that 

p  -  I  dA  cosQ  -  i  cos& 
rz  d.A  rx 

(Ref  3:94),  which  expresses  the  cosine  and  inverse  square  laws  of 
irradiance.  The  inverse  square  law  holds  for  a  point  source  (Ref  3:94). 

We  can  now  derive  an  expression  for  the  reflected  irradiance  of  a  diffuse¬ 
ly  reflecting  flat  plate  as  a  function  of  viewing  and  illumination  angles. 

The  power  incident  upon  a  flat  plate  satellite  at  one  astronomical 
unit  (A.U.),  the  mean  distance  from  the  earth  to  the  sun,  is 

•t^E.Acosp  (Watts) 

where  EQ  is  the  solar  irradiance  in  the  visible  bandpass  at  1  A.U. 
w /  m£),  A  is  plate  area  (m^),  and  ^  is  the  angle  between  the  plate 


normal  and  the  direction  of  the  sun  (Ref  7:15).  The  power  reflected 
from  the  plate  in  the  direction  of  the  sensor  is 

f>  cold  where  f>  is  the  diffuse  reflectiv¬ 

ity  in  the  visible  bandpass,  and  ot  is  the  angle  between  the  plate 
normal  and  the  sensor  line  of  sight  (Ref  7:15).  Combining  these  gives 
the  power  reflected  in  the  direction  of  the  sensor, 

<$f- E.  f  A  Ci>Sp  LOS  at 

The  radiant  exitance  from  the  plate  is  (Ref  3:92), 

M  =  -  =■  CCS  il  CCS  ac 

3  A 

and  radiance,  L  is  given  by 


assuming  a  Lambertian  source.  Since  the  satellite  is  a  point  source 
(Ref  3:536),  reflected  irradiance,  Ep  is 


where  a  is  the  radius  of  a  circular  flat  plate  and  r  is  the  slant  range 
from  the  sensor  (Ref  3:535-536).  Therefore,  the  point  source  reflected 
irradiance  is 


The  portion  of  the  irradiance  equation  which  contains  the  sun-sat- 
ellite-sensor  geometry  is  called  the  "phase  function"  of  the  flat  plate, 
Fp(<*,p).  In  general, 

£t-  F(*,fS,  4>) 

where  ^>is  the  phase  angle.  For  a  flat  plate, 

FP(*,fO=  cosP 

irr* 

This  is  the  simplest  of  the  phase  functions  used  to  determine  the 
integrated  reflected  irradiance  of  the  shapes  used  to  model  the 
diffuse  reflection  characteristics  of  satellites.  Diffuse  irradiance 
equations  for  simple  shapes  are  summarized  in  Table  I 1-2  (Ref  7). 
Sun-Satel 1 i te-Sensor  Geometry 

Phase  Angle,  and  Sensor  and  Solar  Aspect  Angles  In  order  to  use 
the  irradiance  equations  given  in  Table  II-2,  the  sensor  and  solar 
aspect  angles,  otand  |S  respectively,  and  the  phase  angle, f  ,  must  be 
known.  To  determine  and  ^  at  some  instant  in  time,  we  must  know 
the  position  of  the  satellite  in  its  orbit,  the  position  of  the  sensor, 
the  position  of  the  sun  and  the  orientation  of  the  satellites  principal 
axes  within  a  single  coordinate  frame.  Angles  oc,^  and  ^  are  illus¬ 
trated  in  Figure  II-2,  and  Figure  I 1-3  shows  the  three  coordinate  frames 
which  enter  into  the  orbit  prediction  problem. 


To  emphasize  the  effect  that  the  phase  angle  can  have  on  a  signa¬ 
ture,  Figure  I 1-4  depicts  two  satellites  in  high-inclination  low  earth 


XKRAOIANCE 


£  -  Ep/>  ft  COSot  ftfjp 

itr* 


d*  diameter 

Ps=  specular  reflectivity 

/^diffuse  reflectivity 


r  -  Eopdh  sin  m  sinfi[sm6+fH-8)fosdl 
4^*  rL 

~.r  «*> -cost,  cos f  ' 

w  co  L  sm  w  nw  p  J 


4  =  diameter 

h  -  WeiyWt 
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orbits,  visible  simultaneously  to  a  sensor,  S.  The  sensor  Is  In 
darkness  about  one  hour  before  dawn.  The  sensor's  local  horizon 
Intersects  the  plane  of  orbit  2  at  points  C  and  E.  The  arc  OF  marks 
the  Intersection  of  the  plane  of  orbit  2  with  the  cone  of  the  earth's 
shadow.  Satellite  1  (SAT  1 )  Is  sunlit  for  the  entire  time  It  Is  above 
the  sensor  horizon,  but  It  Is  Illuminated  from  behind  and  ^  Is  large 
so  the  observed  diffuse  Irradlance  Is  small.  SAT  2  enters  the  earth's 
shadow  at  point  D,  but  while  It  Is  visible,  ^  Is  small,  so  observed 
diffuse  Irradlance  Is  high,  since  most  of  the  portion  of  the  satellite 
facing  the  sensor  Is  Illuminated. 

Obtaining  Phase  Angle  and  Aspect  Angles  This  paper  uses  the 
universal  variable  formulation  for  time  of  flight  to  solve  the  orbit 
prediction  problem  (Ref  1 : 191-212) -  Radius  and  velocity  vectors,  and 
epoch  time  were  provided  for  each  photometric  signature  received  from 
the  ADIC.  Since  these  vectors  are  epoched  during  the  time  of  track, 
and  the  tracks  are  from  one  to  three  minutes  long.  It  was  not  deemed 
necessary  to  account  for  orbital  perturbations  or  to  advance  the  sun 
position  from  its  position  at  start  of  track.  Sun  positions  were  inter¬ 
polated  from  the  1982  Astronomical  Almanac  (Ref  17).  Astrodynamlc 
constants  are  from  the  DOD  World  Geodetic  System,  1972  (Ref  8).  Table 
I 1-3  lists  applicable  WGS-72  constants. 

The  vectors  obtained  from  the  orbit  prediction  calculation  are  the 
radius  vector,  r  and  the  velocity  vector  9,  in  the  geocentric-inertial 


ASTRODYNAMIC  CONSTANT 


VALUE 


One  Earth  Radius  (Canonical  distance  unit): 
Oblateness  of  the  Reference  Ellipsoid: 
Canonical  Tiro  Init: 

Earth  Gravitational  Parameter: 

Earth  Angular  Rotation  Rate: 


6370.135  km 
.00181881066 
13.44683295  min 
3.986008X10^  km3 


7.292115147X10"4rad 

sec 


Table  I 1-3  WGS-72  Astrodynamic  Constants 
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coordinate  system  (refer  to  Figure  I 1-3) .  We  can  obtain  the  unit  vector 
in  the  direction  of  the  sun  through  a  simple  spherical  rectangular  coor¬ 
dinate  transformation,  if  we  know  its  position  in  right  ascension  (RA) 
and  declination  (DEC).  Right  ascension  is  measured  in  degrees  eastward 
from  the  first  point  of  Aries  along  the  celestial  equator,  and  declina¬ 
tion  is  measured  north  or  south  from  the  celestial  equator. 

s  s  [cOS(otO  cos(P«)Ji+  jeos(06C.)sw(RA)Jj4  £si*(DEOJ  jc 

A  A  A 

where  I,  J  and  K  are  the  unit  vectors  defining  the  axes  of  the  geocentric 
-inertial  coordinate  system. 

The  sensor  position  vector  is  found  using  the  W6S-72  ellipsoidal 
earth  model.  Referring  to  Figure  I 1-5  (Ref  1:95), 


a«  i 

%  =  -7a.il.  ===== — t  h  cosL  . 

J\-t‘  Jl»*t  "  ’ 

g«(t-e)  S1I>L  where 

/i-elswlL 


is  one  equatorial  earth  radius,  e  is  the  eccentricity  of  the  refer¬ 
ence  ellipsoid,  L  Is  sensor  lattitude  and  h  is  the  sensor's  altitude 
above  sea  level  (Ref  1:98).  The  sensor  position  vector,  R  is  given  by 


R  «  (%  tos©)i  +{%  0)j  ♦  ^ 

0 s  0cj  where 
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d  Is  local  sidereal  time,0g  is  Greenwich  sidereal  time,  andA£  is 
the  eastern  longitude  of  the  sensor  (Ref  1:99). 

At  any  point  in  time, 

03=000  ♦(l.0O273'?*093)(nT)(‘D)  radians, 

where  0  is  Greenwich  sidereal  time  at  O*1  U.T.  on  the  first  of  January. 
9o 

D  is  the  number  of  days  which  have  elapsed  since  1  January,  and  1.002737- 
9063  is  one  day  of  mean  solar  time  (Ref  1:101-104). 

The  sensor  line  of  sight  vector  may  be  found  easily,  knowing  r  and 
ft.  From  Figure  II-3, 

i- 

The  remaining  vector  needed  to  determine  ot  and  (3  is  the  one  defin¬ 
ing  the  longitudinal  axis  of  symmetry  of  the  satellite.  Many  payloads 
of  intelligence  interest  have  a  known  nominal  orientation  in  the  body- 
centered  coordinate  frame.  A  coordinate  transformation  from  the  body 
frame  to  the  geocentric- Inertial  frame  provides  the  needed  vector  which 
we  will  call  A=  I  +(C*j)J  . 

To  obtain  the  sensor  aspect  angle  «t,  we  use 


The  solar  aspect  angle  is  given  by 


phase  angle, 4  is  given  by 


4^  180-  cos' 


<fj_j 

VlJII  Si 


Other  Viewing  Geometry  Considerations  The  phase  functions  for 


simple  shapes  given  in  Table  I 1-2  account  for  the  phase  shadowing  of 
the  objects  based  upon  solar  and  sensor  aspect  angles.  When  the  basic 
shapes  are  combined  to  form  a  more  complex  satellite  model,  the  model 
photometric  signature  cannot  be  determined,  in  general,  by  a  simple 
linear  combination  of  the  irradiances  of  components.  The  total  irra- 
diance  is  complicated  by  the  casting  of  shadows  on  some  components  by 
others,  and  by  line-of-sight  obscuration  of  some  parts  by  others. 

These  effects  are  significant  in  diffuse  signatures  when  they  involve 
components  which  have  large  surface  areas,  such  as  solar  paddles. 

Small  protrusions  and  surface  features  are  generally  not  significant 
to  the  diffuse  signature  component. 

This  thesis  models  a  satellite  which  has  two  large,  sun-tracking 
solar  paddles  which  strongly  influence  the  observed  signature.  Figure 
I 1-6  illustrates  a  hypothetical  series  of  satellite  images  as  seen  by 
a  sensor  on  a  pass  for  which  the  paddles  are  illuminated.  As  sensor 


aspect  angle  changes,  the  paddle  away  from  the  sensor  becomes  more 
and  more  obscured  by  the  main  body.  Part  of  the  paddle  closest  to  the 
sensor  also  obscures  part  of  the  cylindrical  body.  The  effects  are 
significant  because  the  surface  areas  involved  are  large. 

The  satellite  model  contained  in  the  subroutine  S1GB1,  described 
in  Chapter  III,  attempts  to  account  for  the  portion  of  the  spaceward 
solar  paddle  which  is  obscured  by  the  main  satellite  body.  The  devel¬ 
opment  of  the  algorithm  for  doing  this  is  given  in  Appendix  B. 


III.  PROGRAM  "SATELLITE  IDENTIFICATION  (SATID)" 
FUNCTIONAL  DESCRIPTION 


Program  SATID  is  written  in  standard  FORTRAN  77  (Ref  6).  The 
coded  listing  appears  in  Appendix  A,  along  with  supplementary  infor¬ 
mation. 

The  Mainline  Program,  SATID: 

Inputs  The  following  quantities  are  required  inputs  to  program 
SATID: 

L. 

1.  Greenwich  sidereal  time  at  0  U.T.  on  1  January  of  the  current 
year.  This  input  is  necessary  only  once  per  year,  and  is  used 
to  calculate  Greenwich  and  local  sidereal  times  at  the  time  of 
track  in  order  to  determine  the  sensor  position  vector. 

2.  The  alphabetic  sensor  code  of  the  sensor  which  collected  the 
data  being  analyzed.  The  codes  used  in  SATID  are  the  ones  in 
current  use  at  the  ADIC.  The  program  uses  sensor  code  to 
select  the  applicable  sensor  longitude,  latitude  and  altitude 
above  sea  level.  Only  the  MOTIF  and  SITU  codes  are  currently 
implemented. 

3.  The  pattern  number  assigned  to  the  data  being  analyzed.  This 
number  is  assigned  to  the  data  at  the  ADIC,  and  permanently 
identifies  the  data  for  later  retrieval  from  permanent  storage. 
SATID  uses  the  number  to  identify  the  data  for  the  analyst  in 
the  output. 


4.  A  geocentric-inertial  radius  vector  in  canonical  units  of  WGS-72 
earth  radii. 

5.  A  velocity  vector  in  WGS-72  earth  radii  per  day.  The  more  usual 
units  of  earth  radii  per  universal  time  unit  were  not  used  because 
vectors  provided  by  the  AOIC  were  in  earth  radii  per  day.  The 
program  converts  to  the  other  units  before  performing  orbit  pre¬ 
dictions. 

6.  The  epoch  time  of  the  input  vectors.  Items  4,5,  and  6  are  used  In 
the  orbit  prediction  calculations. 

7.  The  right  ascension  and  declination  of  the  sun,  in  radians,  at  the 
start  of  track.  These  are  used  to  determine  a  unit  vector  In  the 
direction  of  the  sun  for  Irradlance  calculations. 

8.  The  desired  time  Increment  between  observations,  in  canonical  time 
units.  The  program  can  use  any  time  increment,  but  it  was  kept  at 
one  observation  per  second  for  this  thesis.  The  increment  is  used 
to  update  time  of  flight. 

9.  The  number  of  observations,  or  data  points  to  be  read  Into  the 
program.  This  determines  the  number  of  synthetic  data  points  to 
be  created  by  the  program. 

10.  A  satellite  photometric  signature  in  digital  form. 

Outputs  Program  SAT ID  provides  the  following  outputs: 

1.  Identifying  information  including  the  sensor  by  name,  the  pattern 
number,  and  the  start  time  in  U.T. 


2.  Keplerian  orbital  elements  and  the  right  ascension  and  declination 


of  the  orbit  normal. 

3.  The  following  items  in  parallel  columns: 

a.  Seconds  since  start  of  track 

b.  Azimuth  of  the  satellite  with  respect  to  the  sensor 

c.  Elevation  of  the  satellite  with  respect  to  the  sensor. 

d.  Start  range  of  the  satellite  from  the  sensor  in  kilometers. 
Items  a  through  d  were  used  for  comparison  with  look  angles 
provided  by  NORAD  to  confirm  that  input  radius  and  velocity 
vectors  were  correct  and  that  the  orbit  prediction  algorithms 
worked. 

e.  Right  ascension  of  the  line  of  sight. 

f.  Declination  of  the  line  of  sight.  Items  e  and  f  may  be  used 
to  determine  if  an  unusual  feature  on  the  input  signature  is 
attributable  to  the  satellite  or  to  something  in  the  background 
star  field. 

g.  Phase  angle. 

h.  Four  sets  of  synthetic  data  points  corresponding  to  the  four 
currently  implemented  satellite  models,  SIGA1,  SIGB1,  SIGC1 
and  SIGD1.  Units  are  absolute  visual  magnitudes.  The  term 
absolute  refers  to  normalization  of  all  magnitudes  to  a  range 
of  1000  kilometers. 

i.  A  listing  of  data  points  from  the  true  signature,  TRUSIG. 

4.  A  statistical  results  summary  in  tabular  form.  The  columns  cor¬ 
respond  to  the  four  satellite  models.  The  rows  are,  beginning 

with  row  1: 

a.  The  mean  deviation,  mu,  between  true  and  synthetic  data 
points. 
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b.  The  standard  deviation,  SIGMA,  between  true  and  synthetic 
signatures. 

c.  The  sum  of  the  squares  of  the  residuals,  SSR,  for  each 
satellite  model.  The  model  with  the  minimum  SSR  provides 
the  best  match  to  the  true  data. 

Program  Logic  Figure  III-l  is  a  logic  flow  chart  of  the  mainline 
program  intended  to  clarify  this  narrative  description  of  the  tasks 
performed  by  SATID. 

The  input  satellite  signature  is  loaded  into  a  one-dimensional 
array  of  1000  points  capacity,  TRUSIG.  The  first  major  computational 
task  of  the  program,  after  reading  input  data  and  performing  unit  con¬ 
versions,  is  the  calculation  of  Keplerian  orbital  elements  from  the 
input  vectors.  The  mainline  calls  subroutine  ELSET  which  uses  the 
method  described  in  Bate,  Mueller  and  White  (Ref  1:61-67)  to  compute 
orbital  elements. 

A  large  DO-loop  contains  all  further  computational  tasks  with  the 
exception  of  the  statistical  results  summary  calculations.  The  DO-loop 
control  variable,  N,  is  set  to  the  number  of  observations  in  the  input 
signature,  in  this  case,  the  number  of  seconds  of  track  to  be  analyzed. 
The  loop  counter,  M,  goes  from  1  to  N.  The  following  calculations  are 
performed  for  each  second  of  track  (or  value  of  M)  inside  the  large 
loop: 

1.  Time  of  flight  is  Incremented  from  M-l  to  M  seconds.  If  M*l, 
time  of  flight  is  zero  and  later  computations  are  based  on  the 
original  input  vectors.  The  first  input  data  point  occurs  at 
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delta  time  zero. 


2.  New  radius  and  velocity  vectors  are  calculated  using  the  universal 
variable  formulation  for  time  of  flight  (Ref  1:191-212). 

3.  Current  sensor  position  and  line  of  sight  vectors  are  calculated 
as  described  in  Chapter  II. 

4.  Sensor  look  angles  are  calculated  using  simple  rectangular  spherical 
coordinate  transformations  to  determine  azimuth  and  elevation.  A 
decision  structure  corrects  azimuth  for  the  appropriate  quadrant. 
Range  is  obtained  from  the  magnitude  of  the  line  of  sight  vector. 
Phase  angle  is  the  angle  between  the  previously  determined  line  of 
sight  and  the  sun  vector. 

5.  A  second  DO-loop  within  the  first  contains  the  absolute  visual 
magnitude  calculations  for  each  satellite  model.  The  loop  control 
variable,  Q*  is  set  to  the  number  of  satellite  models  currently 

in  the  program.  In  this  case,  Q'  is  4.  The  loop  counter,  Q,  goes 
from  1  to  Q*.  If  Q  equals  1,  the  program  calls  subroutine  SIGA1. 
SIGA1  stands  for  "signature  of  model  type  Al,"  where  A  is  an 
arbitrarily  assigned  designator  for  a  type  of  soviet  satellite, 
and  1  indicates  it  is  a  "first  generation"  object.  SIGA1  cal¬ 
culates  the  absolute  visual  magnitude  of  the  type  Al  satellite, 
and  stores  the  point  in  an  array  called  SUNA1,  for  "simulated 
signature  of  type  Al."  The  inner  loop  calls  each  satellite 
model  synthetic  signature  generation  subroutine  in  turn,  and 
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a  magnitude  is  calculated  by  each,  and  Is  stored  in  a  corresponding 
simulated  signature  array. 

When  each  satellite  signature  model  subroutine  has  calculated  and 
stored  a  magnitude,  the  flow  returns  to  the  large  DO- loop,  the  loop 
counter,  M  Is  incremented  by  1,  and  the  entire  process  described  in  1 
through  5  above  Is  repeated  for  the  next  second  of  track. 

When  the  large  DO-loop  Is  complete,  the  simulated  signature  arrays 
are  complete  and  the  mainline  program  calls  the  subroutine  COMPAR, 
which  compares  the  true  signature,  point  by  point,  to  each  simulated 
signature  and  outputs  the  statistical  results  summary  table. 

The  Subroutines 

The  following  are  very  brief  subroutine  descriptions  and  the  top- 
level  logic  flow  diagrams  for  some  of  them: 

Subroutine  ELSET  ELSET  stands  for  "element  set,"  and  calculates 
the  Keplerlan  orbital  elements  and  the  right  ascension  and  declination 
of  the  orbit  normal.  The  logic  flow  diagram  Is  Figure  III-2. 

Subroutine  ANGLES  ANGLES  calculates  the  sensor  aspect  angle  (ALPHA) 
and  solar  aspect  angle  (BETA)  of  an  earth-center  stabilized  satellite, 
and  the  sensor  and  solar  aspect  angles  (ALPHAH  and  BETAH)  of  an  horizon- 
stabilized  satellite.  The  earth-center  stabilized  object  has  Its  long 
axis  along  the  orbital  radius  vector,  and  the  horizon-stabilized  object 
has  Its  long  axis  parallel  to  the  vector  formed  by  the  cross  product  of 
the  orbit  normal  and  the  radius  vector.  ANGLES  Is  called  by  the  simu¬ 
lated  signature  generation  subroutines,  for  Irradlance  calculations. 


The  logic  flow  diagram  Is  Figure  III-3. 

Subroutine  ELIPS1  and  2  ELIPS1  calculates  the  Y-coordinates  of 
points  on  the  ellipses  formed  by  the  projection  of  cylinder  endplate 
perimeters  into  the  optical  Image  plane,  corresponding  to  the  X-coor- 
dl nates  of  solar  paddle  corner  points  which  lie  inside  the  sides  of 
cylinder  1  defined  In  Figure  B-l.  ELIPS2  does  the  same  for  cylinder  2. 
These  coordinates  are  used  In  the  plane  analytic  geometry  calculations 
of  subroutine  GE0M1.  The  logic  flow  is  Figure  III-4. 

Subroutine  SE0M1  GE0M1  uses  the  line  of  sight  vector,  sun  vector, 
and  the  orbital  radius  vector,  along  with  type  B1  satellite  dimensions, 
to  define  the  geometry  necessary  to  calculate  the  solar  paddle  area 
which  1$  both  Illuminated  and  visible  to  the  sensor,  even  though  part  of 
one  paddle  Is  obscured  by  the  body  of  the  satellite.  The  approach  Is  to 
establish  an  Image  plane  coordinate  system  with  the  projection  of  the 
radius  vector  defining  the  Y-axis  and  with  the  X-axis  perpendicular  to 
the  Y-axis  and  positive  to  the  sensor's  right.  The  edges  of  all  body 
parts  are  projected  Into  the  Image  plane  and  critical  points,  such  as 
paddle  comers  and  ellipse  centers  are  located.  The  subroutine  also 
locates  the  points  of  Intersection  of  Important  lines  and  ellipses,  and 
creates  equations  for  ellipses,  and  point-slope  form  line  equations. 
GE0M1  Is  called  by  S1GB1.  The  GE0M1  output  Is  used  by  subroutines  CASES 
AREAS,  AREAS1 ,  AREAS2 ,  and  AREAS3. 

Subroutine  Cases  CASES  identifies  which  of  five  viewing  geometries 


applies* to  the  type  B1  satellite,  concerning  the  paddle  obscuration 
calculation.  Case  zero  refers  to  zero  paddle  comer  points  visible  to 
the  sensor.  Case  1  Is  one  comer  point  visible,  and  so  on  to  case  4, 
which  Is  all  four  corners  visible.  CASES  Is  called  by  subroutine 
S1G81.  Figure  111-5  shows  the  logic  flow. 

Subroutines  AREAS1,  2  and  3  These  short  subroutines  evaluate 
Integrals  to  obtain  the  area  between  two  curves.  Limits  of  Integration 
are  defined  In  the  logic  of  subroutine  AREAS,  which  repeatedly  calls 
these  subroutines.  AREAS1  computes  the  area  between  two  lines,  AREAS2, 
between  a  line  and  an  ellipse,  and  AREAS3  between  two  ellipses.  Figure 
III-6  shows  the  logic  flow  for  all  three. 

Subroutine  CP123  CP123  stands  for  "comer  points  one,  two  and 
three."  This  subroutine  performs  an  area  computation  which  must  often 
be  repeated  for  CASE  3,  when  three  corner  points  are  visible.  CP123  Is 
called  by  AREAS.  The  flow  diagram  Is  Figure  III-7. 

Subroutine  AREAS  AREAS  calculates  the  partial  paddle  area  visible 
to  the  sensor  for  a  type  B1  satellite.  The  subroutine  determines  area 
based  upon  the  case  Identified  by  CASES,  and  further  decision  logic 
which  Identifies  subcases  within  each  case.  The  output  Is  the  area  of 
the  paddle,  APAD,  which  Is  actually  the  product  of  the  sensor  aspect 
angle  to  the  paddle  normal  (ALPHAP),  times  the  true  visible  paddle  area. 
AREAS  Is  called  by  subroutine  SIGB1. 

Subroutine  CONE  CONE  approximates  the  diffuse  Irradlance  of  a  cone 


or  a  truncated  cone.  The  conic  is  modeled  by  200  flat  strips  which  are 
triangular  for  a  pure  cone  and  trapezoidal  for  a  truncated  cone.  CONE 
is  called  by  SIGC1,  and  receives  its  inputs  from  the  S1GC1  satellite 
model  parameters,  including  reflectivity,  half-angle,  cone  height,  slant 
length,  base  radius,  and  nose  radius  when  applicable.  The  output  of  CONE 
is  an  approximate  diffuse  Irradiance  for  a  conic  satellite  component. 
Although  S1GC1  is  now  the  only  subroutine  to  call  CONE,  any  future  sub¬ 
routine  could  use  CONE  If  the  appropriate  vector  defining  the  conic  axis 
of  symmetry  is  input  from  the  calling  subroutine.  The  Irradiance  algo¬ 
rithm  is  described  in  Chapter  II.  Figure  III-8  shows  logic  flow. 

Subroutines  SIGA1,  SIGB1,  SIGC1,  and  SIGD1  The  four  satellite  model 
subroutines  calculate  Individual  Irradiances  for  each  component  of  the 
satellite  model  using  the  phase  functions  presented  In  Table  II-2.  The 
component  irradiances  are  summed  and  the  absolute  visual  magnitude  Is 
calculated.  The  magnitude  is  then  read  Into  the  simulated  signature  ar¬ 
ray  for  the  appropriate  satellite  model,  the  array  element  subscript 
corresponding  to  the  current  value  of  the  large  DO-loop  counter  variable, 
M,  described  in  the  paragraph  on  mainline  program  logic. 

All  of  the  satellite  model  subroutines  call  subroutine  ANGLES, 
except  for  SIGD1  which  models  a  sphere  and  requires  only  the  phase  angle, 
which  is  calculated  in  the  mainline.  SIGB1  must  determine  the  observed 
irradiance  of  both  an  unobscured  solar  paddle  and  a  solar  paddle  which 
may  be  partially  obscured  by  the  main  body  of  the  satellite  with  respect 


to  the  sensor.  It  calculates  main  body  irradiance,  unobscured  paddle 
irradiance  and  obscured  paddle  irradiance.  To  obtain  the  latter,  the 
exposed  paddle  area  is  determined  by  subroutine  calls  to  GE0M1,  CASES 
and  AREAS.  SIGC1  contains  a  conic  component,  and  must  therefore  call 
subroutine  CONE.  Logic  flow  diagrams  are  given  in  Figures  III-9 
through  I I 1-12. 

Subroutine  COMPAR  COMPAR  performs  statistical  comparison  of  the 
deviations  between  the  simulated  signatures  and  the  input  true  signa- 
tures.  COMPAR  deals  with  the  simulated  signatures  sequentially, 
beginning  with  SIMA1.  The  subroutine  begins  with  a  large  DO-loop  which 
has  the  number  of  satellite  models,  Q',  as  its  loop  control  variable, 
and  the  model  number,  Q,  as  the  loop  counter.  When  Q=l,  the  program 
compares  the  true  signature  to  SIMA1,  when  Q  equals  2,  the  comparison  is 
with  SIMB1,  and  so  on.  The  simulated  signatures  are  loaded  into  COMPAR's 
working  array,  SIMSIG,  then  each  data  point  in  the  true  signature  array, 
TRUSIG,  is  subtracted  from  its  counterpart  in  SIMSIG. 

dl;  *  Sy.  -  where  s1  is  the  ith  simulated 

data  point,  M^  is  the  ith  measured  data  point,  and  d.  is  the  deviation. 
The  deviations  are  stored  in  an  array  named  DEVSIG.  The  DEVSIG  array 
elements  are  then  summed  and  divided  by  the  number  of  points,  N,  to 
yield  mean  deviation,  mu(/U). 


n 


n 


50 


As  mean  deviations  are  determined  for  each  model,  they  are  stored 


in  array  MEAN.  The  squares  of  the  deviations  in  DEVS1G  are  then  comput 
ed  and  sunned,  and  standard  deviation,  SIGMA  is  calculated,  and  loaded 


Finally,  COMPAR  calculates  the  sum  of  the  squares  of  residuals  for 
each  model,  SSR,  and  loads  them  into  array  SMSR. 


The  arrays  MEAN,  STDEV  and  SMSR  are  strictly  for  output  of  the 
statistical  results  summary  table.  Figure  II 1-13  is  the  COMPAR  logic 
flow  diagram. 
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Figure  I I I- 3  ANGLES  Logic  Flow 
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Figure  I I 1-5-1  CASES  Logic  Flow 


Figure  1 1 1—5—2  CASES  Logic  Flow 


Figure  I I 1-7  CP123  Logic  Flow 
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Figure  111-10  r.IGOl  Logic  Flow 


Figure  1 1 1-12  SIGD1  Logic  FI 


Figure  III- 13-2  COMPAR  Ionic  Flo u 


IV.  SATELLITE  MODELS  AND  VALIDATION  RESULTS 


Model  Descriptions 

Three  types  of  foreign  satellites  were  modeled.  These  are  described 
in  detail  in  the  classified  addendum  to  this  thesis.  To  keep  this  thesis 
unclassified,  none  of  the  actual  signatures  are  identified  by  mission 
class  or  by  the  object  number  assigned  by  the  NORAD  Space  Computational 
Center  (SCC)  in  Colorado  Springs.  The  models  were  given  arbitrary  alpha¬ 
numeric  designators  having  no  relationship  to  similar  designators  employed 
by  the  ADC  Intelligence  Center.  These  designators  appear  within  the 
names  of  the  satellite  model  subroutines.  The  letter  in  the  designator 
refers  to  a  mission  class.  The  number  refers  to  a  variant  of  the  mission 
class.  * 

Model  A1  Model  A1  components  are  a  cylinder  and  a  flat  plate.  The 
cylinder's  longitudinal  axis  remains  aligned  with  the  orbital  radius 
vector,  and  the  components  have  different  reflectivities.  General  con¬ 
figuration  is  given  in  Figure  IV-1A. 

Model  B1  Model  B1  components  include  three  cylinders  and  five  flat 
plates,  two  of  which  represent  sun-tracking  solar  paddles.  Shadowing 
of  components  by  other  components,  and  obscuration  of  one  solar  paddle 
by  the  main  body  are  modeled.  This  is  the  most  complex  of  the  four  mod¬ 
els,  and  the  only  model  validated  using  actual  signature  data.  This 
model  long  axis  also  remains  aligned  with  the  orbital  radius  vector 
(Figure  IV-1B). 


69 


The  purpose  of  model  B1  was  to  develop  and  test  an  algorithm  for 
computing  diffuse  irradiance  of  an  object  with  complex  geometry,  causing 
partial  obscuration  and  shadowing  of  some  parts  by  others. 

Model  Cl  Model  Cl  components  are  two  cylinders,  one  flat  plate  and 
a  cone.  The  body  long  axis  is  aligned  with  the  vector  formed  by  the 
cross  product  of  the  orbital  angular  momentum  vector  and  the  radius 
vector.  If  the  orbit  were  circular,  this  would  be  the  velocity  vector. 
Such  a  configuration  simulates  horizon  stability  in  a  minimum  drag 
configuration.  Figure  IV-1C  gives  general  shape  and  orientation. 

Model  D1  This  model  is  an  arbitrary  sphere  which  does  not  directly 
correspond  to  any  actual  satellite.  The  phase  function  for  this  model 
is  not  restricted  to  diffuse  reflection,  but  includes  the  specular  com¬ 
ponent,  which  is  always  present  for  a  sphere.  The  ratio  of  specular  to 
diffuse  reflectivities  for  this  model  was  chosen  to  be  three  to  one. 
Initial  Model  Configurations 

Estimates  of  physical  and  dynamic  characteristics  of  satellites 
were  provided  by  the  ADC  Intelligence  Center.  All  models  are  simplifi¬ 
cations  of  these  estimates.  From  initial  estimates,  model  B1  was 
refined  by  comparison  with  actual  signatures.  Other  models  are  not 
validated  by  data,  and  represent  only  rough  approximations  to  actual 
configurations. 


Model  B1  Validation 

The  Data  Model  B1  was  validated  using  five  high  quality  photometric 


signatures  collected  by  the  Satellite  Identification  and  Tracking  Unit 
(SITU),  St.  Margaret's,  New  Brunswick,  Canada.  Although  data  from  the 
Maui  Optical  Tracking  and  Identification  facility  (MOTIF)  were  collected, 
the  orbital  radius  and  velocity  vectors  provided  were  invalid,  and  the 
MOTIF  signatures  could  not  be  used. 

The  signatures  provided  by  the  ADIC  were  not  in  digital  form,  but 
were  in  the  form  of  time  versus  stellar  magnitude  plots.  Figure  IV-2 
shows  one  of  the  plots  used  to  validate  model  Bl.  To  make  them  usable 
by  program  SATID,  the  signatures  were  digitized  at  the  ASD  computer 
center,  using  a  Gould  3054  X-Y  recorder  and  the  MODCOMP  classic  data 
analysis  computer.  The  signatures  were  digitized  at  one  point  per 
second  of  time,  or  one  hertz  sampling  rate.  For  near  linear  diffuse 
data,  a  one  to  three  hertz  sampling  rate  is  considered  adequate  (Ref  9: 
15).  The  data  points  were  recorded  on  tape  and  then  output  on  punch 
cards.  Each  punch  card  contained  a  maximum  of  three  data  points,  and 
the  delta  time  from  start  of  track  for  each  point.  Table  IV-1  is  a 
summary  of  signature  data  used  to  validate  model  Bl. 

Procedure  The  validation  procedure  was  to  run  program  SATID  with 
each  of  the  five  signatures,  and  to  alter  model  Bl  input  parameters  be¬ 
tween  runs  in  order  to  minimize  the  SIGB1  SSR  value.  The  reflectivity 
of  the  main  satellite  body  was  the  parameter  changed,  since  dimensions 
were  considered  to  be  fixed  and  solar  paddle  reflectivities  are  well 
known  (Ref  7:18). 


Sources  of  Error  A  nonzero  SSR  for  a  validated  model  is  a  result 
of  many  sources  of  error,  including: 

1.  Inaccurate  model  dimensions,  configuration  or  orientation 

2.  Contributions  to  irradiance  by  unmodeled  features 

3.  Non-uniform  surface  coatings  of  features  which  are  modeled 

4.  The  degree  to  which  true  surfaces  depart  from  the  Lambertian 
assumption. 

5.  Inaccuracies  in  sky  background  corrections  and  sensor  calibrations 

6.  Noise  in  the  data 

7.  Unaccounted  for  X-Y  recorder  bias  introduced  in  the  digitizing 
process 

8.  Differences  among  individual  spacecraft  of  the  same  type,  such 
as  alterations  in  surface  reflectivity  as  the  surface  is  exposed 
to  the  space  environment,  or  differences  when  manufactured. 

The  first  through  the  fourth  sources  of  error  above  could  be 
decreased  somewhat  by  more  elaborate  modeling,  but  since  the  purpose 
of  the  models  is  to  enable  the  software  to  correctly  identify  satellites 
they  need  only  be  sufficiently  accurate  to  allow  an  analyst  to  distin¬ 
guish  between  type:*  with  some  degree  of  confidence. 

Two  categories  of  signatures  emerged  during  the  validation  phase. 
Category  I  signatures  could  be  matched  well  by  the  simulated  signatures 
when  diffuse  reflectivities  were  in  the  range  typical  for  satellite 
materials.  Antireflection  coated  silicon  solar  cells  have  a  diffuse 


reflectivity  of  about  .06,  and  unpolished  aluminum  lies  in  the  .2  to  .3 
range  (Ref  7:18).  The  reflectivities  which  yielded  the  smallest  SSR  for 
model  B1  were  .06  to  .08  for  the  solar  paddles  and  .42  to  .45  for  the 
other  main  body  components.  Specular  reflectivity  for  unpolished  alu- 
minum  is  about  .42  (Ref  7:18),  so  this  value  for  the  diffuse  model  seems 
a  little  high.  The  satellite  may  be  painted,  or  it  may  have  a  slightly 
glossy  surface.  Signatures  which  fall  into  this  category  are  summarized 
in  Table  IV-1.  These  simulated  signatures  exhibit  the  behavior  expected 
using  diffuse  phase  functions,  with  magnitude  dimming  as  phase  angle 
increases.  The  true  signatures  dimmed  with  increasing  phase  angle  also, 
but  at  a  slower  rate,  making  the  synthetic  signatures  too  dim  at  the 
end  of  track.  Tables  IV-2  through  IV-7  are  actual  computer  output  for 
these  signatures. 

The  category  II  signatures  could  not  be  well  matched  by  simulated 
signatures  of  model  Bl.  Reflectivities  required  to  reduce  SSR  were  very 
high,  ranging  from  .6  to  almost  1.0  for  main  body  components.  The  slope 
of  all  simulated  signatures  had  to  follow  the  increasing  phase  angle, 
dimming  magnitude  rule,  but  the  true  data  often  behaved  exactly  the 
opposite,  becoming  slightly  brighter  as  phase  angle  increased.  These 
signatures  are  summarized  in  Table  IV-1.  For  the  sake  of  comparison, 

SSR  values  are  given  for  reflectivities  the  same  as  those  for  category 
I.  Some  category  II  tracks  are  actually  the  last  half  of  long  category 
I  tracks.  Suffix  a  indicates  the  first  half  of  a  track,  and  suffix  b 


the  last  half. 

Table  IV-1  suggests  that  the  critical  factor  in  determining  category 
I  and  category  II  results  is  the  phase  angle.  The  best  matches  of  simu¬ 
lated  and  actual  signatures  occur  for  tracks  with  small  to  moderate  (30 
to  90  degrees)  phase  angles.  The  category  II  tracks  all  have  moderate 
to  large  (90  to  130  degrees)  phase  angles.  A  plausible  explanation  is 
that  the  B1  satellite  is  far  from  being  a  perfect  Lambertian  reflector, 
and  that  its  non-Lambertian  behavior  becomes  more  obvious  as  sensor  as¬ 
pect  angle  diverges  from  normal  incidence.  The  reflective  properties 
of  many  natural  surfaces  are  approximately  Lambertian  near  normal  in¬ 
cidence.  If  the  B1  satellite's  surface  is  glossy  to  some  degree,  we 
would  expect  the  Lambertian  assumption  to  break  down  at  large  phase 
angles. 


sfnsor:  situ*  st  hargarcts*  '.b,  Canada 


DELTA  BACKGROUND  ABSOLUTE  VISUAL  MAGNITUDE 

SEC  A  2  EL  RANGE  HALLS  DECLOS  PHASE  SIGA1  SI6B1  SIGC1  SI6D1  TRUSIG 
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V.  RESULTS  OF  THE  SATELLITE  IDENTIFICATION  EXPERIMENT 


Table  V-l  lists  the  signatures  which  program  SATID  was  used  to 
identify.  These  signatures,  or  portions  of  signatures,  were  not  used  in 
model  validation.  The  table  shows  that  most  results  are  consistent  with 
the  findings  of  the  validation  runs.  Signatures  with  small  to  moderate 
phase  angles  resulted  in  the  smallest  SSR  for  model  Bl,  and  the  slopes 
of  simulated  and  actual  signatures  were  compatible.  Signatures  with 
large  phase  angles  did  not  yield  small  SSR's  for  model  Bl,  but  misi- 
dentified  the  satellite  as  model  Cl. 

The  important  quantity  is  not  the  value  of  the  SSR  for  model  Bl  by 
itself,  since  this  will  vary  with  track  length  and  signature  quality, 
but  the  model  Bl  SSR  value  compared  to  those  of  other  models  for  the 
same  track.  Signature  3845  illustrates  an  important  point  for  any 
photometric  analysis.  Although  the  Bl  SSR  is  quite  low  compared  to 
those  for  A1  and  Dl,  the  Cl  SSR  is  lower  still.  Two  factors  have  prob¬ 
ably  combined  to  cause  the  misidentification.  First,  the  phase  angles 
are  in  the  range  for  which  the  Lambertian  assumption  begins  to  break 
down  significantly  for  the  Bl  satellite.  Second,  the  viewing  and  illu¬ 
mination  geometry  have  created  a  coincidental  situation  in  which  the 
reflectivity-area  product  of  the  diffuse  Cl  model,  and  the  actual  Bl 
satellite  are  almost  the  same.  Tracks  which  undergo  relatively  small 
changes  of  phase  angle,  only  about  ten  degrees  in  this  case,  are  es¬ 
pecially  susceptible  to  such  ambiguity.  The  possibility  of  multiple 


solutions  is  always  present  with  photometric  data.  If  SATID  were  a 
validated,  operational  program,  the  results  for  signature  3845  would  not 
eliminate  model  81  as  the  possible  solution.  Both  low  SSR  models  would 
have  to  be  considered  possibilities,  and  the  ambiguity  would  have  to  be 
resolved  by  using  results  obtained  from  data  collected  later,  other  types 
of  sensors,  or  other  indicators. 
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DELTA  HACKGrOUVO  ABSOLUTE  VISUAL 
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SENSOR:  SITU*  ST  MARGARETS.  \B t  C  Al\ADA 


DELTA  BACKUP  OU'  0  ABSOLUTE  VISUAL  HAGM  TUOE 

SEC  hi  EL  RANGE  RALCS  OtCLOS  PHASE  SIGA1  SIGB1  SIGC1  . [GDI  tpu$: 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 

The  results  of  the  satellite  identification  experiment  show  that  it 
is  feasible  to  perform  pattern  recognition  of  stable  satellites  employ¬ 
ing  diffuse  reflection  models  of  known  satellite  types,  for  geometries 
for  which  the  Lambertian  assumption  is  approximately  true.  Gamache  and 
LaRosa  reached  a  similar  conclusion  (Ref  13). 

Limitations  of  the  Method 

1.  Meaningful  identification  of  satellites  is  achievable  only  when  the 
satellite  being  observed  has  been  modeled.  It  is  possible  for  the  sig¬ 
nature  of  an  unmodeled  object  could  by  chance  exhibit  a  small  SSR  for 
one  of  the  models  on  a  given  pass.  However,  it  is  unlikely  that  an  un¬ 
modeled  object  will  consistently  yield  a  small  SSR  for  a  single  model 
over  several  passes  unless  the  object  is  in  fact  similar  to  the  model. 

2.  The  results  of  the  satellite  identification  experiment  show  that  the 
one  validated  Lambertian  model  is  good  only  for  small  phase  angles. 

Applications  An  operational  program  performing  pattern  recognition 
from  a  model  library  could  be  used  at  sensors  for  early  identification 
and  reporting  to  the  ADIC,  or  at  the  ADIC  itself  by  the  SOI  technician 
on  duty.  The  recognition  process  would  not  depend  heavily  on  the  exper¬ 
ience  of  the  operator.  Multiple  applications  of  a  SATID  type  program 
include: 

1.  Early  mission  identification  of  newly  launched  payloads 

2.  Possible  UCT  identification 


3.  Monitoring  resident  space  objects  for  changes  in 
orientation,  or  configuration,  or  surface  reflec 


tivity  due  to  aging  or  damage. 

4.  Indicating  the  presence  of  a  design  change  in  a 
known  mission  class. 

5.  Determining  that  an  object  which  is  tentatively 
identified  from  other  sources  is  either  not  what 
the  other  sources  say  it  is,  or  that  some  char¬ 
acteristic  is  non-nominal. 

Operational  Recommendations 

1.  Retain  the  PDAM  program  for  use  as  an  aid  to  determining  gross  size 
and  shape  of  unknown  satellites,  but  do  not  try  to  use  it  as  a 
pattern  recognition  tool. 

2.  Develop  another  photometric  analysis  program  for  use  at  the  ADIC, 
with  basic  structure  similar  to  SATID.  Improvements  over  SATID 
should  Include: 

a.  Adding  the  PDAM  preprocessing  module 

b.  Incorporating  a  more  sophisticated  orbit  prediction 
algorithm 

c.  Developing  diffuse  models  for  the  entire  inventory  of 
satellites  of  Interest,  and  validating  them  using  actual 
signatures.  Diffuse  phase  functions  must  not  assume 
perfect  Lambertian  reflection. 


3.  Evaluate  the  utility  of  having  a  SATIO  type  program  available  at 
GEODSS  sensors  as  well  as  the  AOIC. 

Recommendations  for  Future  Research 

1.  An  analysis  is  needed  to  determine  a  correlation  between  the  rela¬ 
tive  magnitudes  of  the  SSR's  of  validated  models  and  correct  iden¬ 
tification  of  satellites  by  program  SATID.  This  will  require  collection 
of  many  more  signatures  of  the  modeled  satellites,  and  validation  of  the 
other  three  models.  Such  an  analysis  would  reveal  how  accurate  the  sat¬ 
ellite  models  must  be  if  correct  identification  can  be  counted  on  for 
some  desirable  percentage  (say  90%)  of  signatures  processed,  assuming  the 
tracks  are  of  modeled  satellites. 

The  small  number  of  signatures  available  to  this  thesis  project  aiade 
such  a  study  impossible,  even  for  model  SIGB1,  but  the  four  satellites 
modeled  were  of  such  different  sizes  and  shapes  that  the  differences  In 
SSR  were  large  in  most  cases  for  which  the  model  was  valid,  and  the  best 
fit  to  data  could  be  easily  selected  by  Inspection.  If  research  In  this 
area  Is  pursued  further,  and  many  more  satellites  are  modeled,  distinct¬ 
ions  between  more  similar  types  may  not  be  as  obvious,  hence  the  need  for 
establishing  quantitative  measures  of  confidence  in  correct  Identification 
for  each  model  Individually. 

2.  Modeling  of  satellites  without  using  the  Lambertian  reflecting  sur¬ 
face  assumptions  should  be  attempted,  so  that  models  remain  valid  at  high 
phase  angles.  An  operational  satellite  Identification  computer  program 
should  retain  its  validity  for  any  viewing  geometry,  so  that  minimal 
detailed  post  analysis  by  the  operator  Is  necessary. 
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Program  SATID  Code  Listing  and 
Supplementary  Material 


VARIABLE  IJAME  LISTS  AND  DEFINITIONS 


The  following  variable  nane  lists  are  arranged  according  to  sub¬ 
routine  in  alphabetical  order.  If  a  variable  occurs  in  more  than  one 
subroutine,  it  is  listed  only  In  the  first  subroutine  In  which  It  Is 
declared.  Refer  to  the  attached  COMMON  Block  nan  to  determine  which 
variables  are  common  to  which  subroutines. 

The  first  column  In  the  lists  qives  the  FORTRAN  variable  name. 

The  next  column  gives  the  variable  type,  Real(R),  Integer(I),  Character 
(CH)  or  Complex(CX).  The  Definition  column  contains  a  short  descrip¬ 
tion  of  the  variable. 
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VARIABLE  NAME  DEFINITIONS 

(SATID) 

NA:n: 

TYPE 

DEFINITION 

ALPHA 

P. 

Sensor  Aspect  Angle(Earth-Center  Stable) 

ALPIIAII 

r 

"  "  "  (Horizon  Stable) 

ALT 

R 

Sensor  Altitude  Above  Sea  Level 

AZ 

R 

Sensor  Line-of-Sight  Azimuth 

I'ETA 

P 

Solar  Aspect  Angle(Earth-Center  Stable) 

rrTAi; 

R 

"  "  "  (Horizon  Stable) 

c 

R 

Series  Expansion  In  Universal  Variable 

Formulation 

non 

I 

Day  Number 

nr  r.  loo 

n 

Declination  of  the  Line-of-Sight 

nccsuti 

R 

Declination  of  the  Sun 

CCLTAT 

R 

Time  Increment  Between  Orbit  Position  Predictions 

HTPX 

P 

Time  Derivative  of  T  WRT  Universal  Variable  X 

LL 

P. 

Sensor  Line-of-Sight  Elevation 

E2 

P 

E-Cornponent  of  Topocentrlc  Radius  Vector 

F 

R 

Function  in  Universal  Variable  Formulation 

FOOT 

R 

Time  Derivative  of  F 

r, 

R 

Function  in  Universal  Variable  Formulation 

COOT 

P 

Time  Derivative  of  G 

r.Hn 

I 

Hour  Part  of  Greenwich  Sidereal  Time  at  0h  UT  on 

1  January  1962 

C«MO 

1 

Minute  Part  of  Above 

GST  R  Greenwich  Sidereal  Time  at  Time  of  Observation 

CETD  R  Greenwich  Sidereal  Time  at  0h  UT  on  1  January 

1002  (Radians) 


NAME 

TYPE 

DEFINITION 

GSO 

R 

Second  Part  of  Above 

IIH 

I 

Hour  Part  of  Track  Time 

I 

R 

I-Component  of  Input  Radius  Vector 

ICOUNT 

I 

Loop  Counter 

I  DOT 

R 

I -Component  of  Input  Velocity  Vector 

I  DOTH 

R 

I-Component  of  Computed  Velocity  Vector 

IDUPDY 

R 

I-Component  of  Input  Velocity  Vector(Earth  radii 
per  Day) 

ISDN 

R 

I-Component  of  Sun  Vector 

ITER 

R 

Iteration  loon  Control  Value 

12 

R 

I-Component  of  Computed  Radius  Vector 

J 

R 

J-Comnonent  of  Input  Radius  Vector 

JDOT 

R 

J-Conponent  of  Input  Velocity  Vector 

JPOT° 

R 

J-Comnonent  of  Computed  Velocity  Vector 

JPLIRDY 

R 

J-Component  of  Input  Velocity  Vector(Carth  radii 
per  Day) 

jsi 

p 

J-Component  of  Sun  Vector 

J2 

n 

J-Component  nf  Computed  Radius  Vector 

i; 

R 

K-Comnonent  of  Input  Radius  Vector 

;niot 

R 

K-Component  of  Input  Velocity  Vector 

>.  'l  ’i 

r* 

K-Comnonont  of  Computed  Velocity  Vector 

1- 

K-Coroponent  of  Input  Velocity  Vector(Earth  radii 
per  Day) 

KCI 1, 

R 

K-Component  of  Sun  Vector 

1 C 

p 

U-Component  of  Computed  Radius  Vector 

LAT 

R 

Sensor  Latitude 

LON 

P 

Sensor  Longitude 

LST 

R 

Local  Sidereal  Time  at  Time  of  Observation 

DEFINITION 


Loop  Counter 

Magnitude  of  Model  1  Signature 

M  M  II  J  11 

II  H  II  ^  U 

N  M  It  ^  II 

Minute  Part  of  Track  Time 
Loop  Control  Variable 
Mean  of  the  Deviations 
Loop  Control  Variable 

II  H  II 

N  II  M 

Phase  Angle 

Signature  Pattern  .‘lumber 

Loop  Counter 

Loop  Control  Variable 

Magnitude  of  the  Orbital  Padlus  Vector 

Right  Ascension  of  the  Llne-of  Sight 

Right  />scens1on  of  the  Sun 

Ratio  of  two  Vector  Components 

Dot  Product  of  Radius  and  Velocity  Vectors 

E-Component  of  Sensor  Position  Vector 

E -Component  of  Llne-of-SIght  Vector 

I -Component  of  Llne-of-SIght  Vector 

J-Component  of  Llne-of-SIght  Vector 

S-Component  of  Llne-of-SIght  Vector 


NAME 


TYPE 


PCFINITIM 


PHOZ 

R 

Z -Component  of  Line-of-Sight  Vector 

RI 

P 

I-Component  of  Sensor  Position  Vector 

RJ 

R 

J-Component  of  Sensor  Position  Vector 

RK 

k 

K-Component  of  Sensor  Position  Vector 

R5 

R 

S-Component  of  the  Sensor  Position  Vector 

PZ 

R 

Z-Component  of  the  Sensor  Position  Vector 

PxO 

R 

Magnitude  of  the  Input  Radius  Vector 

S 

R 

Series  Expansion  in  the  Universal  Variable 
Formulation 

SEC 

R 

Seconds  of  Time 

StiiSOR 

Cll 

Alphabetic  Sensor  Code 

SIGMA 

R 

Standard  Deviation 

SIMA1 

P. 

Model  A1  Signature  Array(500) 

SIMP1 

R 

Model  Rl  Signature  Array(500) 

SIMC1 

R 

Model  Cl  Signature  Array(BOO) 

SIMD1 

R 

Model  D1  Signature  Array(500) 

S!1A 

R 

Semi major  Axis  of  Orbit 

SRNG 

R 

Slant  Range 

SS 

R 

Seconds  Part  of  Track  Time 

SSR 

R 

Sum  of  the  Squares  of  Residuals 

S2 

P. 

S-Component  of  Computed  Radius  Vector 

T 

R 

Time  Plus  Tlme-of-FUght 

TIME 

R 

Tin*  of  Observation 

TOF 

R 

Tlme-of-FUght 

TRUSIG 

R 

True  Signature  Array (1000) 

TSUCti 

R 

Variable  In  Newton  Iteration  for  Universal 
Variable  Determination 
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UATIE 


TYPE 


DEFINITION 


Magnitude  of  Input  Velocity  Vector 

Universal  Variable  for  Tlme-of-Fl 1 ght 

Perpendicular  Distance  Fror>  Earth's  Rotational 
Axis  to  Sensor,  Oblate  Earth  Model 


X  Squared  Divided  by  Semlmajor  Axis  of  Orbit 

Perpendicular  Distance  from  the  Earth's  Equat¬ 
orial  Plane  to  the  Sensor,  Oblate  Earth  Model 

Square  Root  of  Z 

Intermediate  Value  In  Hyperbolic  Orbit  Calculation 
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VARIABLE  NAME  DEFINITIONS 


(ELSET) 


NAME 

TYPE 

DEFINITION 

APGLAT 

R 

ARGUMENT  OF  LATITUDE  AT  EPOCH 

APG PER 

R 

Argument  of  Perigee 

DECON 

R 

Declination  of  the  Orbit  Normal 

ECCEN 

R 

Eccentricity 

EDOTR 

R 

Dot  Product  of  the  Eccentricity  Vector  and  the 
Radius  Vector 

El 

R 

I-Component  of  the  Eccentricity  Vector 

EJ 

R 

J-Component  •  * 

EK 

R 

K-Component  "  M 

HI 

R 

I-Component  of  the  Angular  Momentum  Vector 

HJ 

R 

J-Component  "  "  "  "  " 

UK 

R 

K-Component  "  "  " 

INC 

R 

Orbital  Inclination 

UDOTE 

R 

Dot  Product  of  the  Node  Vector  and  the  Eccen¬ 
tricity  Vector 

HDOTR 

R 

Dot  Product  of  the  Node  Vector  and  the  Radius 
Vector 

HI 

R 

I -Component  of  the  Node  Vector 

IIJ 

R 

J-Component  of  the  Node  Vector  (K-Comp  always  0) 

PER  I  no 

R 

Orbital  Period 

RAAf! 

r 

Right  Ascension  of  the  Ascending  Node 

RAON 

n 

Right  Ascension  of  the  Orbit  Normal 

SLR 

R 

Seml-Latus  Rectum 

riAME 

TYPE 

DEFINITION 

Sf’AXlS 

R 

Semlmajor  Axis  of  Orbit  (km) 

TRANOfl 

R 

True  Anomaly 

TR'JLON 

R 

True  Longitude  at  Epoch 

VC  I 

P 

I -Component  of  Circular  Velocity  Vector 

VCJ 

R 

J-Component  " 

H 

VC  1C 

R 

K-Component  "  "  " 

H 
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VARIABLE  NAME  DEFINITIONS 
(ANGLES) 


na;-;e 

TYPE 

DEFINITION 

rn 

P 

I-Component  of  the  Radius  Vector($ame  as  12 

In  SATID.  Renamed  in  COMMON/VECTR2/) 

LCJ 

P 

Same  as  Above  for  J2 

ECI' 

P 

Same  as  Above  for  K2 

i.nsi 

1! 

I-Component  of  the  Line-of-Siqht  Vector(Same 
as  RliOI  in  SATII).  Renamed  In  COMMON/ VECTR2/) 

I.OSJ 

u 

Same  as  Above  for  RIIOJ 

I.OSIC 

n 

Same  as  Above  for  RliOK 

PHI 

n 

\ 

Phase  AnqIe(Same  as  PHASE  in  SATID.  Renamed 
in  C0WW/VECTR2/) 

Sl-'il 

n 

I-Component  of  the  Sun  Vcctor(Samo  as  ISDN  in 
SATID.  Renamed  In  COMMON/VE CTR2/ ) 

SI.’NJ 

? 

Same  as  Above  for  JSUN 

SUNK 

n 

Same  as  Above  for  KSUN 

VARIABLE  NAME  DEFINITIONS 

(F.LIPS1  AND  CLIPS2) 

nm 

Tvpr 

MLFINITION 

CIC 

R 

Constant  for  Jnner  Edge  Line  Equation (Acronym) 

CLE 

p 

"  "  Lower  "  "  "  ■ 

COE 

R 

"  M  Outer  H  "  M  ■ 

CUE 

R 

•  ”  Uf-per  "  ■  "  * 

HL1 

R 

X-Coord  of  11 Ipse  2  Center 

UL? 

R 

X-Coord  of  Ellipse  3  Center 

MUl 

R 

X-Coord  of  Ellipse  1  Center 

KL1 

R 

Y-Coord  of  Ellipse  2  Center 

KL2 

f! 

Y-Coord  of  Ellipse  3  Center 

Kill 

R 

Y-Coord  of  Ellipse  1  Center 

LENC1 

R 

Length  of  Cylinder  1 

LENC2 

R 

H  N  II  2 

PADLEil 

P 

Paddle  Length 

PADSEP 

R 

Separation  of  Paddle  and  Main  Body 

PADWIP 

R 

Paddle  Width 

PAOC1 

R 

Radius  of  Cylinder  1 

RADC2 

P 

II  N  II  2 

SLOPIO 

R 

Slope  of  the  Jnner  and  Outer  Edge  Equations (Acronym) 

SLOPUL 

P 

"  "  M  Upper  and  Lower  "  M  " 

SMAJ1 

9 

Semlmajor  Axis  of  Ellipses  1  and  2 

SI’AJIS 

R 

Square  of  the  Above 

SI1AJ2 

R 

Semlmajor  Axis  of  Ellipse  3 

s:v\J2r. 

r> 

i ' 

Square  of  the  Above 

fi/UT 

TYPE 

DEFINITION 

71  ' 

SITU 

P 

Semiminor  Axis  of  Ellipses  1  and  2 

-V 

Sf1Itl2 

P 

H  it  ii  m  2 

.  -  „ 

XCP 

R 

X-Coord  of  a  Corner  Point 

i 

XCP1 

P 

X-Coord  of  Corner  Point  1 

XCP1S 

R 

Square  of  the  Above 

XCP  2 

n 

♦  V 

X-Coord  of  Comer  Point  2 

m 

XCP2S 

R 

Square  of  the  Above 

,VCP3 

n 

i\ 

X-Coord  of  Corner  Point  3 

XCP3C 

P. 

Square  of  the  Above 

►  *** 

XCP  4 

R 

X-Coord  of  Comer  Point  4 

XCP4C 

R 

Square  of  the  Above 

1  JU 

XIICE2 

R 

X- Intercept  of  the  Jnner  Edge  with  Ellipse  2 
TAcronym) 

XIIEE3 

P 

Same  as  Above  for  Ellipse  3 

rur.cp 

n 

Same  as  Above  for  Lower  Edge  and  Ellipse  2 

i 

XILCL3 

R 

Same  as  Above  for  Ellipse  3 

XI0EE2 

f> 

Same  as  Above  for  Outer  Edge,  Ellipse  2 

XinLL3 

P 

Same  as  Above  for  Ellipse  3 

XIUEE2 

n 

Same  as  Above  for  Upper  Edge  and  Ellipse  2 

r* 

. 

XII TX  3 

R 

Same  as  AI>ove  for  Ellipse  3 

YCP1 

R 

Y-Coord  of  Comer  Point  1 

YCP2 

R 

N  N  M  N  g 

r# 

- 

YCP3 

P 

N  M  II  M  3 

* .  - 

YCP4 

R 

•  M  M  *  4 

YC2E2 

P 

Y-Coord  of  Intersection  of  Cyl  2  with  Ellipse  2 

r 
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ija:*e 


TYPE 


DEFINITION 


YE1CP 

P 

Y -Coord  on  Ellipse 
of  a  Corner  Point 

1  Corresponding  to  the  X -Coord 

YE1CP1 

p 

II 

II 

H 

II 

N 

for  Comer  Point 

1 

YE1CP2 

P 

II 

N 

H 

II 

N 

N  M 

N 

2 

YE1CP3 

n 

II 

N 

M 

N 

N 

II  H 

M 

3 

YE1CP4 

R 

It 

II 

II 

N 

M 

M  M 

M 

4 

YE2CP1 

R 

II 

H 

II 

•1 

2 

N  M 

N 

1 

YE2CP2 

n 

M 

II 

II 

N 

N 

M  II 

N 

2 

YE2CP3 

R 

N 

M 

M 

M 

N 

M  II 

N 

3 

YE2CP4 

r» 

M 

N 

M 

N 

N 

N  N 

M 

4 

YE3CP 

F: 

II 

II 

N 

N 

3 

•1  M 

H 

YC3CP1 

P 

II 

II 

N 

N 

M 

•1  N 

II 

1 

YF3CP2 

R 

II 

H 

H 

N 

M 

M  II 

N 

o 

YC3CP3 

P 

N 

N 

H 

N 

N 

N  M 

N 

3 

YT3CP4 

p 

N 

II 

H 

II 

M 

N  M 

N 

4 

YIICC1 

R 

Y-Jntercept 

TAcronyro) 

of  the 

Inner  fdge  with  Cylinder  _! 

YIIEC2 

R 

Same  as  Above  for  Cyl  2 

YIIEE2 

R 

Same  as 

Above  for  Ellipse  2 

Y1IEE3 

R 

Same  as 

Above  for  Ellipse  3 

YILEC1 

P 

N 

N 

N 

II 

Lower  Edge  and  Cyl  1 

YILEC2 

R 

N 

H 

H 

N 

N  II 

II 

■  2 

YILEE2 

P 

H 

II 

M 

II 

h  n 

•1 

Ellipse  2 

YILTE3 

n 

II 

N 

N 

II 

N  H 

•1 

•1 

3 

YIDECl 

p 

II 

M 

M 

II 

Outer  " 

H 

Cyl  1 

NAME 


TYPE 


DEFINITION 


YIOEC2 

R 

Y^hitercept  of  the  Outer  Edge  with  Cyl  £ 

YIOCE2 

R 

II  N 

M 

n  n  n 

N 

Ellipse  2 

YIOCE3 

R 

II  N 

II 

M  U  II 

H 

Ellipse  3 

YII'ECl 

R 

II  N 

II 

M  Upper  11 

II 

Cyl  1 

YIIIFC2 

R 

II  II 

M 

it  H  n 

II 

•  2 

Yinrr: 

n 

II  II 

II 

n  n  H 

II 

Ellipse  2 

Yiurr.3 

i? 

II  II 

II 

H  H  H 

II 

"  3 
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VARIAPLE  NAME  DEFINITIONS 
(GEOM1) 


;:a*<e 

TYPE 

DEFINITION 

CP1IC1 

ci: 

Comer  Point 

_1  In  Cylinder  1^  (Acronym) 

CP1IC2 

Cll 

M 

II 

II  II  M  O 

A. 

CP1IC1 

Cl! 

M 

II 

M  "  Ellipse  1 

CP1IE2 

n; 

II 

N 

M  H  M  2 

CP1IL3 

Cll 

N 

H 

II  II  II  2 

CP2IC1 

c:: 

H 

II 

2  "  Cylinder  1 

CP2IC2 

Cl' 

II 

II 

ii  ii  ii  2 

CP2IE1 

ci: 

N 

II 

"  “  Ellipse  1 

CP2IL2 

ci: 

II 

M 

II  II  N  2 

CP? IE 3 

CP 

N 

N 

itiiM  2 

CP3IC1 

Cl! 

II 

N 

3  "  Cylinder  1 

CP3IC2 

Cll 

II 

N 

MU  H  p 

CP3IE1 

Cll 

II 

M 

"  "  Ellipse  1 

CP3II2 

ci; 

II 

H 

ii  ii  u  9 

CP3IL  3 

ci: 

II 

N 

ii  ii  ii  2 

rr-Mci 

Cl; 

II 

N 

4  "  Cylinder  1 

r.p/,icr 

Cl. 

H 

II 

ii  ii  H  2 

rp/:  in 

ci: 

N 

N 

M  "  Ellipse  1 

CP4IL2 

cr 

N 

N 

M  M  M  2 

CP4IE3 

Cl! 

H 

N 

MUM  3 

I1LTAY 

f> 

1*. 

Iteration  Control  Variable 

ETA 

R 

Projected  Solar  Paddle  Tilt  Angle 

.ill 

P 

I-Conponent  of  Normal  Vector  1 

*:iJ 

r> 

l 

J -Component 

of  Ilomal  V.-ctor  1 

:*Af  r 


TYPE 


r-EFiniTiON 


;:1K  r*  K-Component  of  Normal  Victor  1 

.::’I  i.  I-roi.’.;.oimnt  of  formal  Vector  2 

..2,1  !'  J-fonponent  of  .‘Ionia  1  Vector  2 

!!?!'.  K  K-Comoonent  of  Normal  Vector  2 

II3I  2  I-Conponent  of  normal  Vector  3 

l!3J  P  J-Component  of  Normal  Vector  3 

::r,(  "  K-Component  of  Normal  Vector  3 

PAI  P  I-Component  of  Paddle  Axis  Vector 

PA.J  P  J-Component  of  Paddle  Axis  Vector 

PAK  P  K-Component  of  Paddle  Axis  Vector 

PCI  P  I-Component  of  Paddle  Edge  Vector 

PC.]  R  J-Component  of  Paddle  Edge  Vector 

PCK  P  K-Component  of  Saddle  Edge  Vector 

PS I  R  Projected  Solar  Paddle  Rotation  Angle 

XI  P  Angle  Between  Paddle  Edge  and  Line-of-Sight 

XP I '/  R  X-Coord  of  Paddle  Pivot  Point 

YLI.'IE  P  Y-Coord  of  a  Point  on  a  Line 

YPIV  P  Y-Coord  of  Paddle  Pivot  Point 

Angle  Detween  Paddle  Axis  and  Line-of-Sight 


ZETA 


P. 


VARIABLE  HAfT  DEFINITIONS 
(CASES) 


MMC 

TYPE 

DEFINITION 

CASCO 

Cl! 

Case  Zero,  No  Comer  Points  Visible  to 

Sensor 

CA!'C1 

CM 

Case  One,  One  Corner  Point  Visible  to  Sensor 

CASC2 

CM 

Case  Two,  Two  Corner  Points  Visible  to 

Sensor 

CASE  3 

fll 

Case  Three,  Three  Comer  Points  Visible 

to  Sensor 

CASE4 

Cl! 

Case  Four,  All  Comer  Points  Visible  to 

Sensor 

VARIABLE  NAME  DEFINITIONS 


I 

name 

TYPE 

(AREA SI  AND  AREAS2) 

DEFINITIONS 

apcc 

P. 

Area  Between  Two  Ellipse  Curves 

i 

ABLE 

P 

Area  Between  a  Line  and  an  El  1  i 

ARLL 

P. 

Area  Retween  Two  Lines 

CONST 

R 

Constant 

i 

CON ST 1 

Pn 

Constant 

fO.'IST2 

R 

Constant 

K 

B 

Y-Coord  of  an  FI 1  ipse  Center 

i 

K/ 

n 

.  Ellipse  A  Center 

i;r. 

r 

N  ii  ii  H  rj  II 

,-s. 

SLOPE 

R 

Slope  of  a  Line 

SLOPC1 

p 

II  H  M  II 

SLOPE 2 

R 

It  N  H  II 

SKAJA 

R 

Semlmajor  Axis  of  Ellipse  A 

r , 

■ 

s;v\jb 

P. 

N  N  N  H  H 

»*  " 

Sl'IflA 

R 

Semimlnor  "  "  "A 

SI'INR 

R 

N  N  N  H  Q 

* 

:.i 

P 

Limit  of  Integration 

|*7 

X2 

R 

M  N  M 

4 


VARIABLE  NAME  DEFINITIONS 


NAME 

ALPISAP 

.".PAD 


(AREAS) 

TYPE  DEFINITION 

P.  Sensor  Aspect  Angle  for  a  Plate 

P  Area  of  Partially  Otscured  Paddle  Visible  to 

Sensor  Projected  Into  t*ie  Image  Plane 


ARPAD 


p 


Area  of  an  Unobocured  Paddle  Visible  to  Sensor 


VARIAPLE  NAME  DEFINITIONS 


fO'lllIT 


IIAFAMO 


DEFINITION 

Sensor  Aspect  Angle  to  a  Conic  Surface  Element 
Incremental  Surface  Element  Area 
Conic  Base  Radius 

Solar  Aspect  Angle  to  a  Conic  Surface  Element 
I-Component  of  Conic  Normal  Vector 
J-Component  Of  Conic  Normal  Vector 
K-Component  of  Conic  Normal  Vector 
Pure  Cone  Flag 
Cone  Height 
Loop  Counter 
Conic  Half  Angle 
Irradlance  of  the  Cone 
Nose  Radius  of  Truncated  Cone 
Angle  Between  Adjacent  Surface  Element  Normals 
Reflectivity  of  the  Cone 
Slant  Length  of  the  Cone 
Truncated  Cone  Flag 


2 


7ARIAPLF  NAME  DEFINITIONS 
(SIGA1) 


f:AHT 

TYPE 

DEFINITION 

ARLAC 

r 

Cylinder  Projected  Area 

AREAP 

P 

Plate  Projected  Area 

IRCYL 

R 

Cylinder  Irradlance 

IRPLT 

n 

Plate  Irradlance 

LENGTH 

P 

Cylinder  Length 

RADIUS 

R 

Cylinder  Radius 

REFC 

R 

Cylinder  Reflectivity 

REFP 

R 

Plate  Reflectivity 

TI'ETA 

R 

A  Function  of  ALPHA,  BETA  and  PHASE 

I  0) 

I 

1 


143 


VARIABLE  NAME  HFFIN ITIONS 


(SIGP1) 


MM  If 

TYPE 

DEFINITION 

ARE AC 1 

P 

Area  of  Cylinder  1 

AR LAC 2 

R 

II  N  II  O 

L. 

ARCAC3 

n 

ii  ii  H  2 

ARE API 

R 

"  M  Plate  1 

ARLAP2 

P. 

•1  H  HO 

L. 

ARCAP3 

P 

«  M  H  3 

UETDEG 

R 

BETA  In  Degrees 

IRBPDY 

R 

Total  Irradlance  of  Hal n 

IRCYL1 

P 

Cylinder  1  Irradlance 

IP.CYL2 

R 

Cylinder  2  Irradlance 

IRCYL3 

R 

Cylinder  3  Irradlance 

I  I  T  AD 

R 

Irradlance  of  a  Paddle 

IRPLT1 

P 

"  H  Plate  1 

IRPLT2 

R 

N  II  II  2 

IRPI.T3 

n 

i 

ii  ci  ii  2 

RE  TCI 

R 

Reflectivity  of  Cylinder 

REFC2 

R 

N  MM 

RLFC3 

r» 

l\ 

11  N  N 

RFFP1 

R 

"  "  Plate  1 

PEFP? 

R 

"  M  "  2 

Satellite  Body 


1 

2 

3 


3 


3 


VARIABLE  NAME  PEFINITONS 
(SIGC1) 


NAME 

TYPE 

nCFIllITIOH 

IRC1 

n 

i. 

Ir radiance  of  Cylinder  1 

IRC2 

p 

ii  ii  h  2 

IRPLT 

p 

"  *  the  Plate 

LERI 

R 

Length  of  Cylinder  1 

LCH2 

P 

ii  n  "2 

PLTLEN 

R 

M  *  Plate 

PLT'.'IU 

R 

Width  of  Plate 

RADII'S 

P 

Cylinder  Radius 

TlifTA!  I 

P 

A  Function  of  ALPHA!! ,  rFTAII 

VARIALLE  NAME  DEFINITIONS 


NAMC 

rLVrIG 

sr,rir. 

S’TR 

r.Tncv 

mean 

SIIMSQ 


TYPE 


R 


ft 


n 

I  « 


n 

n 


(COMPAR) 

rcriMiTiou 

Array  of  Deviations(r>QO) 

"orking  Array  of  Synthetic  Signature  Data(500) 
Arrav  of  Model  SSR's 
Array  of  ’‘oriel  SICfAA' s(25) 

Array  of  Model  Moan  Deviations 
Sum  of  the  Squares  of  Deviations 
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Appendix  B 

Solar  Paddle  Obscuration  Algorithm  Used  by  Subroutines 
GE0M1,  CASES  and  AREAS 

The  satellite  model  SIGB1  includes  two  sun-tracking  solar  paddles. 
Unless  the  satellite  is  observed  on  or  near  the  meridian  which  passes 
through  the  sensor's  south  point,  one  solar  paddle  will  be  totally  vis¬ 
ible  to  the  sensor,  but  the  other  will  be  partly  blocked  from  view  by 
the  satellite's  main  body.  This  algorithm  first  determines  the  true 
orientation  of  the  satellite  and  solar  paddles  with  respect  to  the 
sensor,  in  three  dimensions  and  then  projects  the  points,  lines  and 
curves  which  define  the  outline  of  the  satellite  into  a  plane  perpen¬ 
dicular  to  the  line  of  sight  vector, J,  referred  to  as  the  optical  image 
plane.  A  two-dimensional  x-y  coordinate  system  is  defined  in  the  optical 
image  plane  with  the  y-axis  always  lying  along  the  projection  into  the 
plane  of  the  orbital  radius  vector,  r,  and  the  x-axis  orthogonal  to  y 
and  positive  to  the  sensor's  right.  The  origin  is  located  at  the  mid¬ 
point  of  the  satellite's  main  body  cylinder.  Figure  B-l  will  be  a  help¬ 
ful  aid  to  visualization  throughout  the  remainder  of  this  algorithm 
description. 

The  projected  exposed  area  calculation  is  performed  after  the  pro¬ 
gram  has  determined  the  x-y  coordinates  of  the  paddle  corners,  equations 
for  the  lines  and  ellipses  which  define  the  satellite  body,  and  the  angles 
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between  the  paddle  edges  and  the  x  and  y-axes.  The  algorithm  proceeds 
as  follows: 

Determine  Key  Vectors 


By  the  time  6E0M1  is  called,  the  line  of  sight,  J,  the  orbital 
radius  vector  r,  and  the  sun  vector  I0R  are  known.  We  must  determine 
the  paddle  axis  vector  PA,  and  the  paddle  edge  vector  Pf: 

PA  is  orthogonal  to  5,  since  s  is  normal  to  the  solar  paddle,  and 
to  r,  since  r  is  the  body  longitudinal  axis  of  symmetry.  Therefore: 


=  sow-  rOl  -  ( n  sun  *  -  son ,  rk  )  J + ( r^sun  j  -sun;  £ 

+  CPA1O& 


PE  is  orthogonal  to  s  and  to  PA,  therefore: 


PE  *  soS  - 


=  (son-  PAk-  PA,  sunOi- (sun- PA  k-  p^son  *)5  4<so*LpAi  -PAi  swn^t 
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Determine  Key  Angles 

Known  angles  at  this  point  are  sensor  aspect  angle, oc  and  solar 
aspect  angle, p  .  We  must  now  find  the  angle  between  the  paddle  axis 
and  the  line  of  sight, 5  .  the  angle  between  the  paddle  edge  vector  and 
the  line  of  sight, £  ,  the  projected  paddle  tilt  angle, v|  ,  and  the  pro¬ 
jected  paddle  rotation  angle, 


S  =■  cos 


■  I  M  PAi-vj,  Pft.-t/xPdi.] 

L  HI  (PHI  J 


£  at  COS 


i;PEi 

♦4.  PEj 

►  4k  PEk  1 

x\ 

n 

j 

C  (zeta)  andg(xi)  will  be  used  to  obtain  paddle  dimensions  pro¬ 
jected  into  the  optical  image  plane.  Obtaining  the  angles  l^(psi)  and 
)|  (eta)  is  less  straightforward. 

' I '*  may  be  viewed  as  the  angle  between  two  planes,  both  containing jf. 
The  first  plane  contains  X  and  r,  and  the  second  plane  contains J  and  TK. 
Since  the  sensor  is  looking  down J,  it  sees  both  planes  edge-on.  The 
lines  of  intersection  of  these  two  planes  with  the  optical  image  plane 
are  the  two  lines  which  we  see  as  the  projected  paddle  axis  and  radius 
vector,  depicted  in  Figure  B-l.  To  obtain^,  the  angle  between  these 


two  vectors  in  the  image  plane,  we  determine  the  angle  between  the  two 
planes  by  finding  their  normals  and  the  angle  between  them. 


N|  =  Z*F  » 


ay 


=  (ii’Tft  -rjioi-cAnt-nJiOA+Min-fiiiifc- 

=  (N.0x  t-(N.j)3  f(NI|,')t- 

Note  that  Rj  is  parallel  to  the  x-axis.  In  like  manner, 

H*  *  f*  *1  =  (Ml L>i  +(Nlj)  J  +  (N20^ 


Therefore,  ^  is  given  by 


<(>  z  COS 


-»  CHU)(N2j)4(Nl))(N2j)  +  ( 

N  |Hz| 


+  (Nik)(N2lQ 

&l 


Y|  may  also  be  viewed  as  the  angle  between  two  planes  containing J. 
The  first  plane  contains?  and  W  and  the  second  contains  J  and  r.  We 
find  the  normal  to  the  plane  containing?  and  PE  thusly: 

Tb=  I*PL  *  (M3i)i  +(N3}>j+(N3|0K. 
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Determine  x-y  Coordinates  of  Paddle  Corners 

First,  locate  the  paddle  pivot  point  depicted  in  Figure  B-l.  The 
x-coordinate  will  be  called  XPIV,  and  the  y-coordinate  will  be  called 
YPIV.  Corner  points  will  have  coordinates  XCP1,  YCP1  through  XCP4, 

YCP4.  These  are  also  the  FORTRAN  variable  names  for  these  points. 

XPIV*  tm£(ri  *p)  where  rl  is  the  main  cylinder 
radius  and  p  is  the  separation  between  the  cylinder  side  and  the  inner 
edge  of  the  paddle.  Note  that  sin  (rl  +p)  is  the  length  in  the  optical 
image  plane  of  the  line  segment  from  the  origin  to  the  pivot  point. 

ypiv=  -cos  +  smCtri+p) 

Next,  we  locate  the  paddle  corners 

XCPl=XPIV-»-S»«rjSmS(YV/z)  where  w  is  paddle  width,  and  sing 
(W/2)  is  the  length  in  the  image  plane  of  the  line  segment  from  the 
pivot  to  corner  point  1. 

YCPt  «VPlV*co$rj  sin  |  tW/z)  and 

ttP2rwmcpi*SMfanS(^)  where  Lp  is  paddle  length, 
YCP2*-cosf  sing(Lp)-*>Vcrt,  H 

*tP4=  WMV -  n  sin  g  (w/z) 

YtP4  =rpiv  -cosrj  sing  (w/z) 

XCP3  =  XCP4 -rsm  ^sinC  (Lr} 

VCP3  =  YCP4  -  cos+  Sing  (  Lf) 
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Determine  the  Ellipse  Equations 

Projections  of  the  cylinder  ends  onto  the  image  plane  form  ellipses. 
The  general  equation  for  an  ellipse  is 

b*  where  h  and  k  are  the  x  and  y-coordi nates 

of  the  ellipse  center,  a  is  the  semi-major  axis,  and  b  is  the  semi-minor 
axis. 

GE0M1  recognizes  three  ellipses,  formed  by  the  top  and  bottom  ends 
of  the  main  cylinder  (ellipses  1  and  2),  and  by  the  bottom  of  the  second 
cylinder  (ellipse  3).  GE0M1  labels  the  ellipse  center  coordinates  in 
the  following  manner: 


HU1  *  the  x-coordlnate  of  the  center  of  ellipse  1 
KU1  *  the  y-coordi nate  of  the  center  of  ellipse  1 
HL1  *  the  x-coordinate  of  the  center  of  ellipse  2 
KL1  *  the  y-coordi nate  of  the  center  of  ellipse  2 
HL2  =  the  x-coordinate  of  the  center  of  ellipse  3 
KL2  *  the  y-coordi nate  of  the  center  of  ellipse  3 


From  Figure  B-l,  it  Is  clear  that  the  x-coordinates  of  all  ellipse 

centers  are  zero.  The  y-coordi nates  are  given  by 

KUI*  Sin  ot  (  Ll/2,) 

»aj=-stnot(Li/2) 

KL2a  -Sin  ((Ll/2)  4  L7) 

where  LI  is  the  length  of  cylinder  1 
and  L2  is  the  length  of  cylinder  2. 
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The  general  ellipse  equation  for  all  values  of  h=0  becomes: 

a*  bl 

and  solving  for  y  gives 

'k 


k (a> - **) 


The  positive  sign  gives  a  y-coordinate  on  the  top  half  of  the 
ellipse,  and  the  minus  sign  gives  a  y-coordinate  on  the  bottom  half  of 
the  ellipse. 

Determine  the  Paddle  Perimeter  Line  Equations 

The  general  equation  of  a  line  in  point-slope  form  is  y=mx+c, 
where  m  is  slope  and  c  is  a  characteristic  constant.  The  slope  of  the 
PE  vector  projection  in  the  image  plane  is, 

mft  «  tan  (IT 

GE0M1  labels  mpg  as  "StOPIO"  for  "slope  of  the  inner  and  outer 
edges."  Similarly, 

Wpa  =  tan(fT+  f) 

=  "SLOPUL",  for  "slope  of  the  upper  and  lower  edges." 
The  constants  for  each  edge  of  the  paddle  are  evaluated  using  the 


point  slope  form,  c*y-mx,  and  setting  x  and  y  equal  to  a  known  point  on 
the  applicable  line  segment.  Expressed  in  FORTRAN  variable  names. 


cte  -  YPiv  -CSi-OPio^(xPlv)  =  cons+aift-t  {or  jywer  ectye.- 

cot*WZ-(*ioPlC>)ClLC(>Z)  ,,  <  gu4er  * 

CU£ s  YCP2-  (St-OPUt-Xlt^PS.)  i  M  upper  •» 

CL6-  YCP4-(SLOPUlKxcP4)  h  >•  Tower  m 

Note  that  if  SL0PUL=0,  then  SLOPIO  is  undefined  and  CUE=YCP2  and  CLE* 

YCP4. 

Determine  the  Intersection  Points  of  Paddle  Edges  and  Cylinder  Sides 
The  lines  through  each  pair  of  corner  points  intersect  the  lines  through 
the  cylinder  sides  at  some  points.  These  points  are  determined  by  using 
the  point-slope  form  of  the  above  line  equations  to  find  y-values  for  x 
equal  to  cylinder  radii.  This  calculation  is  performed  whether  or  not 
the  paddle  edge  actually  crosses  a  cylinder  side. 

Determine  Other  Critical  Points 

Other  points  critical  to  later  area  calculations  are,  1)  y-values 
of  points  on  ellipses  corresponding  to  the  x-coordinates  of  corner  points, 
when  their  absolute  value  Is  less  than  that  of  a  cylinder  radius,  and  2) 
the  Intersections  of  paddle  edge  lines  and  ellipses  under  certain  condi¬ 
tions.  These  points  are  used  to  establish  the  Intervals  of  Integration 
for  area  calculations  In  subroutine  AREAS. 

The  former  set  of  points  Is  found  by  setting  x  equal  to  the  x-coor- 
dinate  of  the  corner  point  and  solving  for  y  in  the  equation  for  the 
ellipse.  For  example.  In  Figure  B-l,  the  y-coordinate  on  the  lower  part 
of  ellipse  2  for  X=XCP4  Is 

where  b^osot(rl),  and 
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a=fl.  Subroutines  ELIPS1  and  ELIPS2  perform  this  calculation.  The 
y-coordinate  in  this  case  is  labelled  YE2CP4  for  "^-coordinate  on 
ellipse  2,  for  corner  £oint  4." 

The  latter  set  of  points  is  calculated  iteratively,  if  it  is 
determined  that  a  corner  point  actually  lies  inside  one  of  the  ellipses. 
For  example,  if 

(ycp4)i<ri2J  VtP4  <  and  VtP4*Ye*cP4 

then  we  know  that  corner  point  4  is  inside  the  lower 
half  of  ellipse  2.  The  iteration  scheme  begins  in  GE0M1  at  statement 
label  75. 

GE0M1  assigns  a  value  of  yes  (Y)  or  no  (N)  to  a  set  of  character 
variables  which  indicate  whether  or  not  corner  points  lie  within  cylind¬ 
er  1,  cylinder  2,  or  any  of  the  three  ellipses,  i.e.,  if  character 
variable  CP4IE2='Y',  "corner  £oint  4  is  jnside  ellipse  2."  See  the 
variable  listings  in  Appendix  A  for  definition  of  this  and  similar 
variable  names  which  are  designed  as  acronyms.  These  character  variables 
are  used  by  subroutines  CASES  and  AREAS  to  determine  which  of  five  geo¬ 
metrical  conditions  (zero  through  4  comer  points  visible  to  the  sensor) 
applies  at  a  point  in  time.  The  case  governs  which  area  calculation 
algorithm  is  used  by  subroutine  AREAS. 


Determine  the  Applicable  Geometrical  Case 


Subroutine  CASES  contains  a  decision  structure  which  assigns  a 
value  of  Yes  (Y)  or  no  (N)  to  character  variables  CASEO,  CASE1,  CASE2, 
CASE3,  and  CASE  4.  The  number  In  each  variable  name  refers  to  the 
number  of  paddle  corner  points  visible  to  the  sensor. 

CASEO  means  all  corner  points  lie  inside  a  cylinder  or  an  ellipse 
in  the  optical  image  plane,  and  none  are  visible  to  the  sensor. 

CASE1  has  two  subcases.  These  either  indicate  that  corner  point 
two  is  visible  and  all  the  others  are  not;  or  that  corner  point  three 
is  visible  and  all  the  others  are  not. 

CASE2  covers  three  subcases.  Corners  1  and  2  can  be  visible  and 
the  other  not,  or  corners  2  and  3  can  be  visible  and  the  others  not, 
or  corners  3  and  4  can  be  visible  and  the  others  not. 

CASE3  has  two  subcases.  Only  corner  point  one  can  be  Inside  a 
cylinder  or  an  ellipse  and  the  others  visible,  or  only  corners  point  4 
can  be  insider  or  cylinder  or  an  ellipse  and  the  others  visible. 

CASE4  indicates  that  all  corner  points  are  visible. 

Determine  the  Projected  Paddle  Surface  and  Visible  to  the  Sensor 


A  sample  calculation  for  the  geometry  in  Figure  B-l  will  best  illus 
trate  the  method,  Figure  B-l  is  an  example  of  CASE3.  To  determine  the 
exposed  paddle  area,  it  Is  necessary  to  Integrate  to  obtain  the  area 
enclose  by  exposed  paddle  boundaries.  The  area  between  two  functions 
f(x)  and  g(x),  over  an  Interval  (t1tt»)  is  given  by 


Figure  B-l  requires  Integration  over  four  intervals.  Interval  1 
Is  from  tj=r2  to  t2arl,  f(x)  is  the  equation  of  the  line  through  comer 
points  3  and  4,  and  g(x)  is  the  equation  of  ellipse  2.  Therefore, 


Interval  2  is  from  tj=rl  to  t2aXCPl,  f(x)  is  the  equation  of  the 
line  through  corner  points  1  and  4  and  g(x)  is  the  equation  of  the  line 
through  comer  points  3  and  4,  so  similarly,  interval  3  is  from  tj«XCPl 
to  t2=XCP3,  so  that 


Interval  3  is  from  t^XCPl  to  t2=XCP3. 
i  l*tP3 

A3-  ffK*  *  'X(cue)  -  wpq*2 

2  I  Z 


Therefore, 

"XCcLt) 


Finally,  interval  4  Is  from  XCP3  to  XCP2,  giving  us 

-  A(eoe) 


A4  -  2l  4  %( cue) 
2 


MP3 


Total  exposed  projected  paddle  area  is  then 


MP3 

xtPl 

pm 

MP3 


Ap  -  At  *v  +■  A3  "t-  A4 
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Note  that  since  this  is  a  projected  area,  it  would  not  have  to  be 
multiplied  by  cosoc  in  the  irradiance  equation,  and  irradiance  is  given 
by 


where  r  is  slant  range 


The  decision  structure  of  subroutine  AREAS  determines  the  proper 
interval  of  integration  for  each  case  and  subcase,  based  on  the  output 
of  subroutines  GE0M1  and  CASES,  and  calls  subroutines  AREAS1,  AREAS2, 
and  AREAS3  to  calculate  projected  paddle  areas.  The  area  determined 
is  named  APAD,  and  is  used  by  subroutine  SIGB1  in  the  calculation  of 
total  irradiance  of  the  satellite. 


Appendix  C 

The  Subroutine  CONE  Conic  Irradiance  Approximation  Algorithm 


Because  of  the  extreme  complexity  of  the  programming  decision 
structure  that  would  be  necessary  to  evaluate  the  true  conic  phase 
function  (Ref  7:  ),  an  iterative  approximation  is  used  by  subroutine 

CONE  to  obtain  diffuse  conic  irradiance.  CONE  will  approximate  irra¬ 
diance  for  a  pure  cone  or  for  a  truncated  cone. 

The  conic  surface  is  approximated  by  200  flat  strips  as  depicted  in 
Figure  C-l.  In  Figure  C-l,  YC  is  the  circular  velocity  vector,  £N,  is 
the  first  normal  vector  to  a  flat  surface  element  which  is  calculated, 
and  EN.  is  an  arbitrary  ith  surface  element  normal.  The  anglesatj  and 
|S.  are  tiie  sensor  and  solar  aspect  angles  measured  from  the  normal  to 
the  ith  surface  element.  The  angle  is  the  angle  between  adjacent  sur 
face  element  normals.  For  the  number  of  surface  elements, n  =200, <J=1. 8 
degrees.  The  angle  Y  is  the  conic  half  angle.  The  circular  velocity, 
VC,  calculated  by  subroutine  ANGLES,  is  given  by 
VC  -  H  XT 

where  H  is  the  orbital  angular  momentum  vector.  The  body-centered  coor 

A  A  A 

dinate  system  is  defined  by  the  unit  vectors,  r,  vc  and  H.  From  Figure 

A  A  A 

C-2  we  see  that  unit  vector  CN^,  which  lies  in  the  r,  vc  plane,  is 
given  by 


cn,  ~  (cx>s}0r  +(&\nfc)vc  +  (o)h 


We  can  express  CNt  in  terms  of  the  geocentric-inertial  frame,  since  we 
already  know  r  and  vc  in  that  frame. 

?=(ri)X+(rj)J  +(rw)!c  and 

VC  =(VCi)i+(VCj)J  +  (VCy^iL. 

Therefore, 

04,=  cosX(ri+fivffc)+  smfcfVQ+VCj+VCfe.)  and 

A 

geocentric-inertial  components  of  CN|  are  given  by 

CN,i=a>sX(fiH  s\r\\ f(vct) 

CN,j  =  cosX(rj)+smif(YCj)  t  and 

CNih  -  cos*  (r„) + s \f&  (VC,,) 

A 

CN,  is  the  only  one  of  the  flat  surface  element  rrtwmals  wvUh  is 
in  the- r,  vc  plane  at  an  angle  If  from  r.  To  compute  the  components  of 
the  other  199  normal  vectors,  we  perform  a  right-handed  rotation  of  the 
body-frame  about  vc  by  the  incremental  angled),  to  establish  another 
vector,  r',  which  will  be  an  angle  If  from  the  next  normal  vector,  CN2. 
CN2  will  lie  in  the  r  ,  vc  plane,  r '  is  given  by 

A 

r'~  co5u>r-  smtoH 

Geocentric-inertial  components  of  r  are, 

r{ = cos  tAiCn)-  s»n  w C Hi) 
r/=  cosu>(  fj)-  smujfrtj) 
r '  =  cos  -  smu)(Hk) 


Finally,  CNg  is  given  by 


CNi  =  (td£  $  )  f  #  and 

CN2-=  cos*  [c^>su>(riHs^nu)(Hi^5m^f(vCi,) 
CN2jrcdsK  |cosu>(fjHiMu>(Hj)]  +  si*  WvCj) 
cn^-cos*  ["cosu)Cr^-FS\yu*)(H^Jt^«nXfvCn.) 


In  subroutine  CONE,  the  calculation  of  new  normal  vectors  is  done 
in  a  DO-loop  which  runs  from  1  to  200.  Each  time  a  normal  vector  is 


A  /  A 

determined,  <j  is  added  to  the  previous  ti)  and  r  and  CN..  are  calculated 
until  a  normal  has  been  determined  for  each  flat  surface  element.  At 
each  step,  the  sensor  and  solar  aspect  angles  are  calculated.  If  both 
or  and  pi  are  less  than  90  degrees,  an  irradiance  is  calculated  for  the 
surface  element  and  it  is  added  to  total  conic  irradiance.  If  either 
o<j  or  p.  is  greater  than  90  degrees,  no  irradiance  is  calculated  for 
that  step  in  the  loop.  Total  irradiance  is  given  by 


200 


where  E..  is  the  irradiance  of  the  ith  surface 


element  and  Ec  is  total  conic  irradiance,  called  IRCONE. 

Subroutine  CONE  is  currently  configured  for  a  conic  which  is  axi¬ 
ally  aligned  along  vc,  but  it  could  easily  be  made  general  by  allowing 


it  to  receive  a  different  vector  as  input. 
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